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AUTONOMOUS  INTEGRATED  GPS/ INS  NAVIGATION  FILTER  FOR  ADVANCED 

SPACECRAFT  APPLICATIONS 


Triveni  N.  Upadhyay,  George  J.  Priovolos,  Harley  Rhodehamel 
MAYFLOWER  COMMUNICATIONS  CO.,  Reading,  MA 

and 

A.  Wayne  Deaton 

NASA  MARSHALL  SPACE  FLIGHT  CENTER 


The  paper  summarizes  the  results  of  a  recently  completed 
feasibility  study  on  GPS-based  techniques  for  navigation  and 
attitude  update.  The  motivation  for  autonomous  navigation, 
and  some  results  for  both  low-earth  and  high  altitude  users 
are  provided.  The  need  for  autonomy  will  grow  in  the  future 
as  missions  require  faster  decision  making,  continuous 
coverage  and  increased  survivability.  The  GPS/INS 
integrated  navigation  system  offers  one  solution  with 
increased  accuracy  and  mission  flexibility  compared  to 
conventional  systems. 

Three  GPS-based  attitude  determination  techniques  were 
analyzed  for  providing  improved  total  navigation  solution. 
The  autonomous  integrated  navigation  filter  will  use  a 
combination  of  these  techniques,  depending  on  the  spacecraft 
mission  phase,  to  provide  the  best  solution  and  flexibility 
in  mission  planning.  The  reconfigurability  of  the 
integrated  navigation  filter,  in  real-time,  to  adapt  to  the 
spacecraft  mission  phases  is  also  discussed.  A  knowledge- 
based  navigation  filter  manager  (i.e.,  expert  system)  is 
discussed  which  will  be  designed  to  automatically  select  the 
navigation  filter  states  and  filter  modes  depending  on  the 
mission  phase  and  available  resources,  e.g.,  thrust  versus 
coast  phase,  one  GPS  antenna  versus  multiple  GPS  antennae, 
translation  versus  rotational  dynamics,  gyro  failure  versus 
accelerometer  failure. 

The  integrated  navigation  filter  results  presented  in  this 
paper  correspond  to  a  17-state  Kalman  filter  which  was 
implemented  in  its  U-D  factor  formulation  for  improved 
numerical  stability.  Two  spacecraft  applications:  (1)  low- 
earth  orbit  users  (e.g.,  OMV,  Space  Station,  Shuttle-C)  and 
(2)  high  altitude  users  (e.g.,  NASA  GEO  platform,  Lunar- 
return  and  Mars-return  transfer  vehicles)  are  analyzed  and 
numerical  simulation  results  are  presented  for  these  cases. 
The  LEO  user  simulated  for  performance  evaluation  consisted 
of  an  OMV  high-thrust  trajectory  burn  of  approximately  300 
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two  antennas  and  to  IMUs.  We  looked  at  two  cases,  one  where 
GPS  was  available  prior  to  a  burn  phase  (good  initial 
conditions)  and  one  where  GPS  was  not  available  (poor 
initial  conditions) .  Details  are  presented  for  the  second 
case. 

The  simulation  results  indicate  the  excellent  performance  of 
the  filter  even  for  the  case  of  poor  initial  conditions. 

This  test  case  assumed  that  the  OMV  did  not  use  GPS  during 
the  coast  phase  and  a  large  position  and  velocity  error  had 
accumulated  prior  to  the  start  of  the  burn.  The  application 
of  the  integrated  GPS/INS  navigation  filter  allowed  the 
large  initial  condition  errors  to  be  reduced  in  a  short-time 
and  resulted  in  excellent  position,  velocity  and  attitude 
estimation  solution  accuracy. 

The  application  of  the  integrated  navigation  filter  for  a 
high  altitude  user  consisted  of  a  spacecraft  at  GEO 
altitude.  Since  the  GPS  availability  is  limited  at 
altitudes  higher  than  10,000  nmi  (the  GPS  satellite  antennae 
are  pointed  towards  the  earth) ,  the  availability  of  a  high 
gain  antenna  and  a  rubidium  clock  onboard  the  spacecraft  was 
assumed  for  this  analysis.  The  simulation  results  for  a  GEO 
STV  user,  with  the  above  assumptions,  were  obtained.  These 
results  further  support  the  utility  of  GPS  for  high  altitude 
spacecraft  missions.  It  is  believed  that  for  these  missions 
GPS  may  be  the  sole  provider  of  the  required  navigation 
accuracies. 

The  application  of  autonomous  integrated  GPS/INS  system 
technology  is  shown  in  this  paper  to  provide  excellent  total 
navigation  performance,  flexibility  in  contingency  mission 
planning,  and  may  help  eliminate  attitude  update  sensors 
which  will  reduce  the  cost  and  size,  power  and  weight  for  a 
wide  class  of  spacecraft  missions. 

A  successful  application  of  GPS  tracking  of  advanced 
spacecraft  will  alleviate  the  requirement  for  ground 
tracking  (or  TDRSS  tracking),  and  the  real-time  onboard 
reconfigurability  feature  of  the  integrated  navigation 
filter  will  reduce  the  cost  of  mission  planning  and  will 
support  the  NASA  requirement  of  developing  advanced  fault- 
tolerant  spacecraft  missions. 
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AUTONOMOUS  INTEGRATED  GPS/INS  NAVIGATION 
PROCESSING  FOR  THE  OMV 

(Proof-Of-Concept  Demonstration) 
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r\  „  ,  LOW  GAIN  ANTENNA  ASSEMBLY 

6MV  (REVISED  PLUME  IMPINGEMENT  CONFIGURATION 
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Excellent  performance  of  filter  under  both  conditions 


OMV  Simulation  Results 
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EVALUATION  FOR  A  12  HOUR  PERIOD  -  ZENITH  ANGLE  =  110  DEG 
INCLINATI0N=55  DEG,  ALTITUDE=1000  N.MI. 
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AUTONOMOUS  INTEGRATED  GPS/INS  NAVIGATION  FILTER 


I  L  T  (DEGREES)  NORTH  TILT  (DEGREES)  ERST  TILT  (DEGREES) 
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AUTONOMOUS  INTEGRATED  GPS/INS  NAVIGATION  FILTER 
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AUTONOMOUS  INTEGRATED  GPS/INS  NAVIGATION  FILTER 
GPS  Satellite  Visibility  for  the  STV 
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limited  GPS  availability 

STV  (geosynchronous  user  with  almost  poor  GPS  satellite  visibility) 

—  RSS  Position  error:  47  m  (1  -cr) 

—  RSS  Velocity  error:  0.013  m/sec  (1-cr) 
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AN  EXPERIMENT  IN  ATTITUDE,  POSITION  AND  VELOCITY  DETERMINATION 

WITH  ROGUE  GPS  RECEIVERS 

Tom  K.  Meehan 
JET  PROPULSION  LABORATORY 


An  experiment  was  performed  on  an  aircraft  to  investigate  the  attitude  deter¬ 
mination  capability  of  GPS.  This  was  done  in  anticipation  of  an  application 
on  EOS,  to  be  launched  in  1996,  which  will  have  three  antennas  as  shown  and  a 
GPS  receiver  similar  to  the  JPL  Rogue  receiver  with  1  cm  accuracy  for  pseudo¬ 
range  and  submillimeter  for  carrier  phase.  The  goal  with  EOS  is  to  get  real¬ 
time  attitude  determination  of  about  0.1  mrad,  real-time  position  error  of 
about  1  m,  and  post-processing  errors  of  0.03  mrad  and  a  few  centimeters, 
respectively. 

Another  objective  of  the  aircraft  experiment  was  to  prove  the  concept  of  using 
Rogue  GPS  receivers  to  calibrate  SAR  (Synthetic  Aperture  Radar)  measurements 
to  obtain  ocean  current  velocities.  Two  SAR  antennas  on  the  bottom  of  the 
aircraft  fuselage  are  side  looking  (45°),  one  transmits,  both  receive.  From 
the  Interferometric  recombination  of  the  two  received  signals  one  can  deter¬ 
mined  the  ocean  wave  propagation  (i.e.,  current  velocity).  One  significant 
error  source  Is  yaw  and  pitch  in  the  aircraft;  e.g.,  an  error  of  0.01°  in  yaw 
yields  an  error  of  4  cm/s  In  current  velocity  versus  a  requirement  of  about 
0.6  cm/s. 

GPS  was  used  to  demonstrate  the  potential  for  controlling  these  errors  -  the 
initial  goal  was  to  achieve  an  accuracy  of  0.1 0  In  pitch  and  yaw.  Two  anten¬ 
nas  were  mounted  on  a  DC-8  and  two  Rogue  GPS  receivers  were  used  (8  satel¬ 
lites,  LI,  L2  carrier  phases,  PI,  P2  codes,  1  s  data  rate).  A  twenty-minute 
segment  In  which  4  satellites  were  simultaneously  visible  by  both  antennas  was 
extracted  from  the  test  flight.  The  first  graph  shows  the  heading  of  the 
aircraft  as  determined  from  the  velocity  which  was  obtained  by  differencing 
GPS  position  Information.  The  second  graph  shows  the  difference  between  head¬ 
ing  and  azimuth  and  reveals  the  3°  bias  due  to  the  second  GPS  antenna  not 
being  on  the  center  line  of  the  aircraft  and  a  drift  due  to  change  in  wind 
velocity  (pilot  had  to  Introduce  yaw).  Most  of  the  detailed  features  of  this 
plot  are  real  (not  noise)  and  due  to  aircraft  maneuvers. 

The  results  were  better  than  our  goals  for  pitch  and  yaw.  The  largest  error 
source  precluding  the  achievement  of  SAR  requirements  Is  multipath.  A 
repositioning  of  the  rear  antenna  would  Improve  this  as  would  the  availability 
of  more  GPS  satellites. 
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GPS  Antenna  Placement  for 
Eos-A  Platform 
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Mid-Ocean  Mapping  Using  SAR 
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Azimuth]  -  (Aircraft  Hooding]  (Ooqt 


FTCURE  1 


AZIMUTH  Of  AIRCRAFT  VELOCITY 
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Future  Developments 
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PRELIMINARY  GPS  POINTING  DATA  RESULTS 


BY 

PHIL  WARD 

SENIOR  MEMBER  TECHNICAL  STAFF 
TEXAS  INSTRUMENTS  INCORPORATED 
6600  CHASE  OAKS.  BOULEVARD 
P.O.  BOX  869305  MS  8449 
PLANO,  TEXAS  75086 
(214)  575-3824 
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PRELIMINARY  GPS  POINTING  DATA  RESULTS 
ABSTRACT 


GPS  attitude  measurement  promises  to  become  a  cost  effective 
source  of  precise  alignment  and  calibration  for  inertial  systems 
on  space  platforms  whose  orbits  are  within  GPS  coverage  such  as 
low  earth  orbit  satellites,  the  Space  Shuttle  and  the  Space 
Station.  GPS  has  already  proven  itself  to  be  effective  for 
determining  an  ultra  precise  vector  between  two  GPS  antennas  for 
first  order  geodetic  surveying  using  carrier  doppler  phase 
I nterferometry  techniques.  The  same  principles  can  be  applied  to 
obtain  two  orthogonal  pointing  vectors  to  provide  attitude 
determination.  Two  advantages  of  GPS  pointing  over  geodetic 
surveying  are  that  the  antenna  phase  center  separation  (the 
absolute  value  of  the  unknown  vector)  is  known  a  priori  and  the 
antennas  are  one  to  ten  meters  apart  rather  than  thousands  of 
meters.  The  known  antenna  separation  can  be  used  to  speed  up  the 
solution  of  the  phase  ambiguity  problem  while  the  close  proximity 
of  the  antennas  ensures  negligible  differential  atmospheric  delay 
error,  permits  the  differential  multipath  error  to  be  smaller  and 
a  sinqle  dual  port  GPS  receiver  can  be  used  to  obtain  both  sets  of 
observab 1 es . 

This  paper  describes  some  preliminary  field  test  pointing  data 
taken  with  a  Texas  Instruments  Incorporated  (TI)  GPS  Receiver 
which  has  been  slightly  modified  for  this  experiment.  The 
receiver  is  the  new  AN/PSN-9  5-Channel,  C/A-Code,  16  pound, 
portable  manpack/vehicular  (MP/V)  receiver  now  under  production 
contract  to  the  U.S.  Air  Force  Joint  Program  Office.  At  the  heart 
of  the  AN/PSN-9  is  a  receiver  channe 1 -on-a-ch i p  design  which 
embodies  a  TI  Invention  specifically  designed  to  obtain  pointing 
observables  from  two  GPS  antennas,  but  requiring  only  one  GPS 
receiver:  a  5-Channel,  two  port  per  channel  design.  A  tripod 
mounted  beam  with  one  meter  antenna  separation  was  used  to  obtain 
pointing  data  from  real  GPS  satellites  while  continuing  to 
navigate  in  the  normal  manner  with  respect  to  the  reference 
antenna.  TI  has  also  Invented  a  technique  which  resolves  the 
phase  ambiguity  problem  without  external  aiding  within  a  few 
seconds  under  dynamic  conditions  rather  than  the  half  hour  or  more 
typically  required  for  stationary  geodetic  observations.  A  rapid 
convergence  example  is  presented.  Single  difference  and  double 
difference  raw  data  from  real  GPS  satellites  are  illustrated  as 
well  as  the  pointing  results  determined  from  these  observations 
(azimuth  and  elevation)  which  verify  that  the  invention  actually 
works.  The  initial  field  data  taken  with  1  meter  antenna 
separation  indicate  that  an  azimuth  accuracy  of  0.5  milliradlan 
will  be  obtained  with  I  meter  antenna  separation  with  no  special 
antenna  phase  matching,  calibration  or  smoothing  of  the  raw  data. 
Owing  to  the  nature  of  the  TI  invention,  the  noise  contribution 
within  the  receiver  is  essentially  thermal  noise.  The  major 
sources  of  pointing  error  are  expected  to  be  differential 
multipath  and  differential  antenna  phase  center  migration,  but 
further  experimentation  will  be  required  to  characterize  all 
sources  of  error  and  to  determine  what  is  the  limit  to  which  these 
sources  of  error  can  be  reduced.  A  "truth"  pointing  reference  is 
currently  under  construction  and  phase  matched  antennas  will  be 
used  in  future  tests.  These  results  will  be  reported  in  a  second 
paper  at  the  IEEE  PLANS  '90  in  Las  Veqas,  Nevada,  March  20  to  23, 
1990. 
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Artist' 8  concept  of  TI  420  GPU 


T!  420  POSITION/NAVIGATION  SYSTEMS 
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WHAT  IS  UNIQUE  ABOUT  THE  Tl  POINTING  INVENTION? 
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GEODETIC  QUALITY  RECEIVER/OSCILLATOR  NOT  REQUIRED 


PRELIMINARY  GPS  POINTING  DATA 
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T1420  GPU  -  PARKING  LO 

1  SEC  DIFF  PHASE  -  SV  1 2-  7/21/89 
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GPS  TIME  IN  SECONDS 
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TI420  GPU  -  PARKING  LOT  7/21/89 
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DOUBLE  DIFFERENCE  PSD.  METERS**2/HZ 
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PHASE  :  SV  9/1  1 


TI420  GPU  -  PARKING  LOT  7/21/89 


SPECTRUM:  PHASE  SV  9/1 1 
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USE  PHASE  MATCHED  GPS  ANTENNAS 
MINIMIZE  DIFFERENTIAL  MULTIPATH 
ALIGN  AN  INERTIAL  ATTITUDE  PLATFORM 


Tl  GPS  POINTING  INVENTION  WORKS  AS  PREDICTED 


co 

h- 

z 
111 
2 
a : 

LU 

X 

X 

LU 

< 

< 

o 

X 

H 

3 

DC 

h- 

h- 

O 

3 

Q 

Z 

O 

o 

o 

h* 

Z 

3 

0. 


o 

111 

0. 

X 

LU 

Q 

Z 

< 

>~ 

o 

< 

X 

3 

o 

o 

< 

0 


z 

o 

CL 


o 

CO 

LU 

< 

X 

H- 

CO 


z 

O 

X 

< 

X 

LU 

CO 

< 

z 

z 

LU 

I- 

z 

< 

X 

LU 

h~ 

LU 


a 

x 

3 

O 

s 

X 

H 

3 


g  9 

Bi 

h-  in 

a 6 


o 

o 

X 

LU 

Q 


h 

CO 


LU 

Q 


X 

o 

X 

< 

2 


CO 

< 

o 

LU 

> 

CO 


b  5 


o 

05 

CO 

z 


LU 

LU 

LU 


5 

3 

CO 

LU 

X 

h 

X 

o 

X 

UJ 

X 


0 

z 

H 

O 

3 

Q 

O 

O 

CO 


LU 

O 

X 

X 

O 

LU 

> 


II 

g  LU 

co  m 
2  u- 

m  s 

“  O 
CO  o 

b  P 
o  < 
x  E 
q  b 

<;  CO 

O  111 
H  x 

is 

O  rf) 


H  2 

g  X 

n  m 
o  x 

S  W 

X  G 


LU 

Q 

X 

> 

< 

X 


2 

CO 

X 

0 


302 


ALGORITHMS  FOR  SPACECRAFT  ATTITUDE  DETERMINATION  WITH  GPS 

Duncan  B.  Cox  and  Haywood  S.  Satz 
MAYFLOWER  COMMUNICATIONS  COMPANY,  INC. 

and 

Ronald  L.  Beard  and  G,  Paul  Landis 
NAVAL  RESEARCH  LABORATORY 


A  filter  is  under  development  that  utilizes  GPS  data  and  accelerometer  and 
gyro  data,  if  available,  and/or  spacecraft  models  to  estimate  position,  veloc¬ 
ity,  and  attitude,  among  other  variables.  The  emphasis  here  is  on  using  GPS 
data,  in  particular,  the  accumulated  carrier  phase,  not  only  for  attitude 
determination  but  also  for  navigation.  The  difference  between  the  accumulated 
phase  data  and  the  delta  phase  measurements  is  that  in  many  conventional  fil¬ 
ters  that  treat  delta  phase,  the  time  intervals  spanning  the  delta  phase  do 
not  cover  the  entire  time  history  and  data  are  lost.  The  filter  design  has  an 
additional  goal  of  fault  detection  and  Isolation. 

The  Integrated  navigation  filter,  designed  for  autonomy  (real-time  solutions), 
is  based  on  the  Kalman  filter  which  processes  accumulated  phase,  pseudorange, 
delta  phase,  delta  attitude,  and  delta  velocity.  In  the  accumulated  phase 
model,  the  initial  unknown  constant  Is  treated  as  a  state  in  the  filter  (this 
yields  many  states  representing  measurements  of  different  satellites  at  dif¬ 
ferent  antennas) . 

The  simulation  results  are  based  on  observations  to  four  GPS  satellites  using 
three  antennas  (12  simultaneous  measurements)  with  a  6  s  data  rate.  Only 
accumulated  phase  measurements  are  used  and  Initial  position  and  clock  errors 
are  assumed  to  be  15  m;  while  Initial  velocity  is  assumed  to  be  known  to  1  m/s 
and  orientation  error  is  0.1  rad.  In  this  case  the  vehicle  is  not  rotated  and 
ambiguities  (Individual  and  differential)  are  not  known.  Note,  however,  that 
even  position  is  being  estimated  since  Information  Is  coming  from  the  changing 
gravity  gradient  as  the  satellite  orbits  in  space.  Velocity  is  determined 
after  two  samples  because  the  constants  are  assumed  to  be  really  constant  and 
the  accumulated  phase  at  two  points  yields  the  velocity  directly.  Attitude  Is 
not  determined  at  all  because  there  Is  no  rotation.  Rotation  of  the  space¬ 
craft  supplies  Information  on  the  ambiguity  constants  and  attitude  can  be 
determined.  Knowledge  of  the  ambiguity  constants  now  provides  consistency 
checks  of  the  data  (detection  of  cycle  slips). 

In  conclusion,  the  general  concept  for  an  Integrated  navigation  filter  has 
been  developed  with  the  unique  additional  feature  of  accommodating  accumulated 
phase  measurements.  Using  the  accumulated  phase  to  determine  attitude  and 
perform  navigation  is  believed  to  be  superior  to  using  delta  phase.  However, 
more  study  and  simulation  analyses  are  required. 
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SIX-DEGREE-OF  FREEDOM  NAVIGATION  ALGORITHMS 
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GINSS  OVERVIEW 

GPS-Inertial  Sj^stem  Simulation 
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PRELIMINARY  RESULTS  FOR  NON-ROTATING  VEHICLE 

WITH  FREE  AMBIGUITY  STATES 
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PRELIMINARY  RESULTS  FOR  ROTATING  VEHICLE 
WITH  FREE  AMBIGUITY  STATES 
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ROTATING  SPACECRAFT,  UNCONSTRAINED  AMBIGUITIES 


MAYFLOWER  COMMUNICATIONS  COMPANY 


PRELIMINARY  RESULTS  FOR  NON-ROTATING  VEHICLE 
WITH  CONSTRAINED  AMBIGUITY  STATES 
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NONROTATING  SPACECRAFT,  CONSTRAINED  AMBIGUITIES 


Cr 

Ci 


98T»ld  J3  S  «  d«A  S££  Hid  itl  i  «  BttA  Seg“lild  Jt)  £  *  UbA 


320 


MAYFLOWER  COMMUNICATIONS  COMPANY,  INC. 


ALGORITHMS  FOR  SPACECRAFT  ATTITUDE 
DETERMINATION  WITH  GPS 
CONCL  US  IONS 
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Preliminary  Experimental  Performance  of  the  TOPEX  Global  Positioning 
System  Demonstration  Receiver  (GPSDR) 

Lance  Carson 

MOTOROLA  GOVERNMENT  ELECTRONICS  GROUP 

This  paper  summarizes  the  basic  requirements  of  the  TOPEX  Demonstration  Recei¬ 
ver,  overviews  its  architecture,  reviews  the  measurement  techniques,  and  of¬ 
fers  some  preliminary  results.  The  TOPEX  receiver  is  a  fully  space-qualified 
receiver  designed  and  built  from  the  ground  up.  The  key  requirements  imposed 
by  JPL  for  the  TOPEX  receiver  include  continuity  in  measurements  of  LI  and  L2 
carrier  phase  simultaneously  on  6  channels;  autonomous  selection  of  GPS  satel¬ 
lites  that  are  optimum  for  TOPEX;  high  accuracy  in  time,  carrier  phase,  and 
pseudorange,  and  low  interchannel  bias;  and  the  ability  to  use  an  existing 
external  frequency  reference.  The  architecture  of  the  receiver  implements  a 
double  conversion  for  LI  and  L2  with  the  local  oscillator  as  common  reference 
both  for  maximum  symmetry  and  delay  tracking,  and  samples  are  digitized  at  40 
MHz  rate.  Channel -on-a-chip  technology  is  used  (6  chips  for  each  frequency)  - 
one  difference  in  the  TOPEX  chip  architecture  is  that  measurements  are  inte¬ 
grated  to  the  20  msec  level  allowing  the  processor  to  accommodate  many  other 
functions,  such  as  navigation  and  interface  software. 

The  program  status  shows  that  hardware  and  software  design  are  completed  and 
that  the  engineering  model  hardware  has  been  built.  Hardware/software  inte¬ 
gration  is  in  progress  and  the  two  flight  receivers  for  TOPEX  are  ready  to  be 
assembled. 

The  technique  used  to  measure  the  random  error  of  carrier  phase  and  pseudo¬ 
range  is  based  on  the  difference  of  two  channels  that  are  tracking  two  differ¬ 
ent  satellites  having  the  same  dynamics  (this  can  be  done  with  a  simulator), 
thus  eliminating  drift  and  oscillator  errors.  The  performance  of  the  advanced 
development  model  receiver  is  close  to  the  predicted  value  of  about  1  mm  for 
the  carrier  phase  and  about  14.7  cm  (versus  8  cm  predicted)  for  the  pseudo¬ 
range.  The  latter  disparity  may  be  due  to  some  residual  contribution  from  the 
simulator  entering  even  after  the  differencing  technique  is  applied.  By 
tracking  just  one  satellite  with  all  channels  yields  the  interchannel  biases 
which  turn  out  to  be  less  than  the  noise  of  the  measurements. 

Discussion: 

Question:  Are  there  any  data  on  the  reliability  (l.e.,  MTBF)  of  the  TOPEX 
receiver? 

Answer:  We  have  some  rough  numbers  based  on  correlations  with  similar  equip¬ 
ment  (on  the  order  of  100000  hours). 

Question:  Are  the  Interchannel  biases  constants? 

Answer:  They  seem  to  be  more  random  in  this  particular  case;  but  Ideally  with 
the  digital  architecture  of  the  receiver  they  should  be  zero.  We  are  still 
investigating  the  mechanism  for  these  biases. 

Question:  Are  the  chips  dedicated  to  process  either  LI  or  L2,  or  both? 

Answer:  They  are  Identical  but  their  physical  location  contralns  them  to 
process  either  LI  or  L2. 
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TOPEX  GPSDR  KEY  REQUIREMENTS  SUMMARY 
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DESIGN  AND  PERFORMANCE  OF  THE 
GPS  RECEIVER  UNIT  (GPSRU)  FOR  THE 
NASA  ORBITAL  MANEUVERING  VEHICLE 


Roger  M.  Weninger 
Richard  Sfeir 

Satellite  &  Space  Electronics  Division 
Rockwell  International 
Anaheim,  California 


ABSTRACT 

This  paper  provides  a  description  of  the 
spaceborne  GPS  Receiver  Unit  (GPSRU) 
planned  for  installation  on  NASA's 
Orbital  Maneuvering  Vehicle  (OMV).  The 
receiver's  hardware  and  software  imple¬ 
mentation,  and  its  interfaces  with  the 
OMV,  will  be  presented.  Simulation 
results  demonstrating  the  receiver's 
performance  will  also  be  reported.  A 
breadboard  unit  has  been  operating  at 
Rockwell  since  May  1989.  The  qualifica¬ 
tion  test  is  planned  for  August  1991, 
and  the  flight  units  are  scheduled  for 
delivery  in  March  1992. 


INTRODUCTION 

In  August  1988  a  contract  was  awarded  to 
the  Satellite  and  Space  Electronics 
Division  of  Rockwell  International  to 
provide  GPS  receivers  for  the  Orbital 
Maneuvering  Vehicle  (OMV).  This  recei¬ 
ver  is  based  on  a  Miniature  GPS  Receiver 
(MGR)  chipset  developed  by  Rockwell 
under  contract  to  DARPA  (1).  This  chip 
set  is  shown  in  Figure  1  and  their  phys¬ 
ical  parameters  are  summarized  in 
Table  1.  It  consists  of  a  GaAs  MMIC 
RF/IF  Translator,  a  Silicon  Frequency 
Synthesizer,  a  CMOS  Signal  Processor,  a 
CMOS  Multi-Function  Interface  and  an 


Figure  1.  Miniature  GPS  Receiver  (MGR)  Chip  Set 
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Table  1.  MGR  Chip  Set  Summary 


CHIP  TYPE 

f - 

DEVICE 

COUNT 

DIMENSIONS 

(inches) 

IMPLEMENTATION 

TECHNOLOGY 

POWER 

(mU) 

RF/IF  Translator 

420 

0.340  x  0.350 

Gallium  Arsenide 

1500 

Signal  Processor 

19,800 

0.185  x  0.220 

1.25  Micron  Bulk  CMOS 

90 

Multifunction  Interface 
(MFI) 

29,000 

0.370  x  0.370 

1.6  Micron  Bulk  CMOS 

20 

Advanced  Architecture 
Microprocessor  (AAMP) 

74,000 

0.250  x  0.261 

2  Micron  Bulk  CMOS 

90 

Frequency  Synthesizer 

350 

0.165  x  0.165 

Bipolar  Silicon 

500 

Advanced  Architecture  Micro-Processor.  b)  Use  the  LI  signal  of  1.57542  GHz 

This  chip  set  is  configured  with  memory,  only, 

a  frequency  standard,  a  power  supply, 

and  user  interface  circuitry  to  make  up  c)  Input  signal  power  of  the  LI  C/A 

the  OMV  GPSRU  described  herein.  The  signal  within  the  range  of  -138  dBm 

receiver  will  have  two  channels  capable  to  -104  dBm. 

of  tracking  either  C/A  or  P  code,  and 

will  determine  position,  velocity  and  d)  Achieve  navigation  accuracy  within 

time  for  the  OMV.  Two  receivers  will  be  10  minutes  (90%  probability)  after 

installed  on  the  spacecraft,  each  an  initialization  accuracy  of: 

operates  independently  of  the  other. 

time  +  1  second 

The  OMV  prime  contractor  is  TRW  and  the  position  +  80  nmi 

responsible  NASA  agency  is  the  Marshall  velocity  +  600  ft/sec 

Space  Flight  Center. 

e)  Pseudo  range  error  not  to  exceed  28 

SPECIFICATION  REQUIREMENTS/  feet  for  C/A-code  and  6  feet  for 

PREDICTED  PERFORMANCE  P-Code  (1  sigma). 

The  following  performance  parameters  f)  Pseudo  range-rate  error  not  to 

were  extracted  from  the  OMV  GPSRU  spe-  exceed  0.08  feet/second  (1  sigma), 

cification  document  (2).  Actual 

receiver  performance  is  expected  to  be  g)  Absolute  GMT  time  error  of  1  micro¬ 
better  than  these  values.  As  an  example  second  (1  sigma), 

the  simulated  navigation  position  and 

velocity  errors  summarized  in  Table  3  h)  Navigation  position  error  per  axis 

are  better  than  the  specified  per-  not  to  exceed  397  feet  (3  sigma) 

formance  defined  in  h)  and  i).  when  GDOP  a  4.6. 

a)  Minimum  carrier  to  noise  ratio  of  i)  Navigation  velocity  error  per  axis 

33.2  dB-Hz  for  the  C/A  signal  and  not  to  exceed  1.19  feet/second  (3 

30.2  dB-Hz  for  the  P-code  signal.  sigma)  when  GDOP  =  4.6. 
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j)  The  receiver  shall  be  capable  of 
correcting  the  corrupted  data 
resulting  from  Selective  Availa¬ 
bility  (SA). 

The  specified  (2)  system  errors  for  the 
evaluation  of  navigation  performance 
are  shown  in  Table  2.  Using  these 
errors  in  a  computer  simulation  to 
predict  GPSRU  navigation  performance 
produced  the  results  shown  in  Table  3. 
These  results  show  a  predicted  one- 
sigma,  spherical  position  error  with 
P-code  of  11.7,  11.9  and  12.6  feet  for 
OMV  altitudes  of  125,  200  and  1000 

nautical  miles,  respectively. 


FUNCTIONAL  DESCRIPTION 

The  two  channel  GPSRU  sequentially 
tracks  a  selected  group  of  GPS  satel¬ 
lites.  The  channels  function  independ¬ 
ently  so  that  data  from  two  satellites 
can  be  processed  in  parallel.  During 
normal  tracking,  the  two  channels  are 


allocated  a  one-second  processing  inter¬ 
val  for  each  satellite,  i.e.,  within  one 
second  all  of  the  real-time  functions 
associated  with  acquisition,  tracking 
and  navigation  are  accomplished.  At  the 
end  of  the  one  second,  the  two  channels 
reacquire  and  track  the  next  two  satel¬ 
lites  in  their  respective  satellite 
assignments . 

The  receiver  can  operate  on  either  or 
both  the  Clear  Access  (C/A)  code,  or  on 
Precision  Code  (P-code)  depending  upon 
performance  requirements.  As  a  general 
rule  when  estimated  position  is  not  well 
known  C/A  code  is  used  for  acquisition 
and  then  handover  occurs  to  the  P-code; 
however,  if  the  estimated  position  and 
time  error  is  within  +300  P-code  chips 
(+30000  feet) ,  a  satellite  can  be 
acquired  with  P-code. 

A  functional  configuration  of  the  GPSRU 
is  shown  in  the  flow  diagram  of  Figure  2 
which  is  designed  to  show  the  signal/ 
data  flow  from  the  antenna  input  to  the 


Table  2.  GPS  System  Error  Assignment 


PSEUDO  RANGE  ERROR 

1-SIGMA  (ft) 

ERROR  SOURCE 

P-CODE 

C/A  CODE 

Space  Segment: 

Clock  and  Navigation  Sub-system  Stability 

9.8 

9.8 

Prediction  of  Satellite  Perturbation 

3.3 

3.3 

Other 

1.6 

1.6 

Control  Segment: 

Ephemerls  Prediction  Model  Implementation 

13.8 

13.8 

Other 

3.0 

3.0 

Tropospheric  Delay  Compensation 

0 

0 

Multipath 

3.9 

3.9 

Other 

1.6 

1.6 

Receiver  Noise  and  Resolution 

6.0 

28.0 

Ionospheric  Delay  Compensation 

17.0 

17.0 

RSS 

25.5 

37.4 
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Table  3.  Predicted  Navigation  Performance 


ALT. 

OF 

OMV 

(nmi) 

NO. 

OF 

GPS 

SATS 

CODE 

TYPE 

lo  POSITION 
ERROR  (ft) 

lo  VELOCITY 

ERROR  (ft/sec) 

PROBABILITY  ERROR 
EXCEEDS  3o  LIMITS 

lo 

CLOCK 

ERROR 

(ns) 

P0STII0N 

VELOCITY 

X 

D 

Z 

X 

Y 

Z 

X 

D 

X 

Y 

Z 

125 

21 

P 

6.5 

7.0 

6.8 

0.040 

0.046 

0.045 

0 

0 

0 

0 

0 

0 

5.0 

125 

21 

C/A 

7.9 

8.6 

9.1 

0.039 

0.043 

0.042 

0 

0 

0 

0 

0 

0 

5.8 

200 

21 

P 

6.3 

7.9 

6.3 

0.039 

0.044 

0.042 

0 

0 

0 

0 

0 

0 

4.9 

200 

21 

C/A 

■ 

8.9 

7.5 

0.039 

0.046 

0.043 

0 

0 

0 

0 

0 

0 

5.6 

1000 

21 

P 

n 

7.5 

6.7 

0.039 

0.045 

0.048 

0 

0 

0 

0 

0 

0 

5.2 

1000 

21 

C/A 

8.5 

9.0 

9.1 

0.044 

0.052 

0.045 

0 

0 

0 

0 

0 

0 

5.5 

data  output  to  the  host  vehicle.  The 
flow  diagram  also  shows  the  functional 
partitioning  of  receiver  elements  among 
RF  and  digital  hardware,  and  software 
processing. 

RF  PROCESSING 

The  incoming  Li  signal  is  routed 
through  an  overload  protection  limiter 
with  broadband  low  signal  linearity  and 
a  low  noise  amplifier  which  is  followed 
by  a  band  limiting  filter.  The  signal 
is  then  down-converted  to  an  IF  of 
173.91  MHz  by  the  injection  of  1401.51 
MHz  signal  to  the  first  mixer.  The  down 
converted  signal  is  further  filtered 
and  amplified  through  an  IF  bandpass 
filter  and  an  IF  amplifier.  The  IF 
signal  is  then  fed  to  an  inphase  and  a 
quadrature  mixer  to  generate  the  I  and 
Q  signals.  The  I  and  Q  mixer  local 
oscillator  frequency  is  obtained  by 
dividing  the  1401.51  MHz  signal  by  8, 
yielding  a  frequency  of  175.18875  MHz. 
The  resultant  I  and  Q  mixer  outputs  are 
video  signals  at  a  frequency  of  1.27875 
MHz.  Each  of  the  I  and  Q  signals  are 
followed  by  a  low  pass  filter  and  a 
wideband  video  amplifier.  The  video 


amplifier  is  gain  controlled  over  a 
range  of  30  dB.  These  stages  are  used 
for  Automatic  Gain  Control  (AGC). 

The  1.27875  MHz  IS.Q  AGC  output  signals 
are  converted  to  digital  signals  in  A/D 
converters.  The  A/D  conversion  process 
is  accomplished  using  a  null  zone  limi¬ 
ter  which  provides  a  +1,  0,  or  -1  output 
as  a  function  of  input  signal  amplitude. 
The  null  zone  limiter  was  chosen  over  a 
binary  hard  limiter  since  it  maintains 
the  signal  processing  simplicity  while 
providing  superior  performance.  It 
should  be  noted  that  the  finite  quanti¬ 
zation  of  the  A  to  D  converter  will  be 
reflected  in  a  processing  loss  rather 
than  a  measurement  error  or  bias.  The  I 
and  Q  quantized  signals  are  fed  to  a 
pair  of  CMOS  drivers  to  drive  the  sig¬ 
nal  processing  functions. 


SIGNAL  PROCESSING 

The  major  functions  performed  by  signal 
processing  are  carrier  demodulation  and 
C/A  or  P-code  removal.  As  illustrated 
in  Figure  2,  each  of  the  two  receiver 
channels  has  its  own  signal  processor 
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Figure  2.  GPSRU  Functional  Flow  Diagram 
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chip.  Each  signal  processor  synchro¬ 
nous1  samples  the  incoming  I  and  Q 
sigt.s  at  twice  the  P-code  rate  (20.46 
MHz;  and  inputs  them  to  the  carrier  and 
code  tracking  function. 

The  sampled  I  and  Q  signal  is  a  1.27875 
MHz  carrier  shifted  by  any  doppler  fre¬ 
quency  and  quadrature  phase  modulated 
by  the  1.023  MHz  C/A  code  modulo-two 
summed  with  the  50  Hz  navigation  data, 
and  in-phase  modulated  by  the  10.23  MHz 
P-code  also  modulo-two  summed  with  the 
50  Hz  data.  The  I  and  Q  samples  are 
input  to  a  digital  phase  rotator  driven 
by  the  carrier  VCO  output.  The  phase 
rotator,  acting  as  a  single  sideband 
mixer,  translates  the  signal  to  base¬ 
band  by  removing  the  carrier  and  its 
doppler  frequency  leaving  the  sampled 
C/A  code,  P-code,  and  data. 

The  signal  is  then  passed  to  code 
removal  logic  which  correlates  the 
selected  code  (C/A  or  P-code)  *  from  the 
Code  Generator  to  the  incoming  code.  The 
Code  Generator  output  is  essentially 
moved  back  and  forth  in  time  until  it 
exactly  matches  the  incoming  code,  at 
which  point  the  code  signals  cancel  and 
only  the  data  plus  any  tracking  errors 
remain.  The  amount  of  time  shift 
required  to  make  the  codes  match  is 
directly  proportional  to  pseudorange. 
An  Integrator  then  sums  the  sampled  data 
signal  into  1.0  millisecond  samples 
coherent  with  the  assumed  C/A  epochs 
for  software  processing. 

CHANNEL  PROCESSING 

Channel  processing  functions  derive 
carrier  phase,  code  delay,  and  signal 
amplitude  estimates  from  the  raw  (1  ms) 
signal  samples.  The  processed  signals 
are  further  summed  and  directed  -  to.  the 
code  and  carrier  loop  filters,  which 
are  integrated  at  a  50  Hz  (20  ms)  rate 
and  provide  control  outputs  to  the 
hardware  VCOs.  Although  primarily  con¬ 
cerned  with  signal  tracking,  channel 
processing  provides  signal  search, 
acquisition,  reacquisition  and  subframe 


data  demodulation  as  well.  Internally 
the  channel  processing  function  accepts 
signal  samples  at  a  1  kHz  rate  from  the 
signal  processor  and  provides  back  com¬ 
mands  at  a  50  Hz  rate.  The  channel 
processing  function  also  provides  the 
pseudo-range  measurement  from  the  code 
loop,  the  delta  pseudo-range  measurement 
from  the  carrier  loop,  and  the  carrier 
power  to  noise  spectral  density  ratio 
(C/N0)  estimate.  It  interfaces  to  the 
navigation  and  interface  processing 
functions  at  a  1  Hz  rate,  providing 
channel  measurements  and  status,  down¬ 
link  subframe  data  (when  required),  and 
accepting  channel  command  data. 

NAVIGATION 

The  navigation  function  accepts  pseudo¬ 
range  and  delta  pseudorange  measure¬ 
ments  ,  combines  them  with  the  corres¬ 
ponding  GPS  satellite  position,  velo¬ 
city,  and  acceleration  in  Earth- 
Centered-  Inertial  (ECI)  coordinates, 
and  calculates  optimal  estimates  of  the 
host  vehicle  position,  velocity  and 
time.  Linearized  measurement  data  and 
measurement  error  variances  are  input 
to  an  8-state,  extended  Kalman  filter 
which  estimates  corrections  to  the 
nominal  navigation  solution.  The  filter 
also  maintains  estimated  error  statis¬ 
tics  in  the  form  of  an  error  covariance 
matrix. 

In  addition  to  providing  navigation 
outputs,  the  navigation  solution  pro¬ 
vides  estimates  of  pseudorange,  pseudo¬ 
range  rate  and  pseudorange  acceleration 
to  the  code  correlator  to  aid  the 
receiver  during  acquisition  to  position 
the  Code  Generator  for  sequential 
tracking  operation. 

INTERFACE 

The  GPSRU  has  an  Instrumentation  Port 
(IP)  that  provides  a  standard  bidirec¬ 
tional  serial  data  interface,  accepts 
data  capture  commands  and  generates  a 
precise  time  mark  at  a  1  Hz  rate.  The 
IP  is  used  for  testing  and  is  suitable 
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as  a  host  vehicle  interface.  The  data 
format  is  compatible  with  ICD-GPS-204 
data  format  for  an  Instrumentation  Port. 
This  is  compatible  with  an  RS232C  serial 
interface  port  operating  at  9600  baud. 

The  GPSRU  also  has  an  OMV  specified 
interface  consisting  of  discretes  and 
of  four  buffers  with  unidirectional 
serial  interfaces.  There  are  three  out¬ 
put  buffers  and  one  input  buffer  to 
provide  for  the  following  data  transfer 
and  receiver  control: 

Input  Buffer: 

Mode  control  (Off,  Standby  1,  Standby 
2,  Built-In-Test  and  Navigate) 

Initialization  data  (position,  velo¬ 
city  and  time)  in  ECI 

GPS  constellation  almanac  data 

Thrust  vector  (when  available)  in  ECI 

Update  output  data  buffer  requests 

Output  Buffer  1: 

X,  Y  and  Z  position  in  ECI  coordinates 

X,  Y  and  Z  velocity  in  ECI  coordinates 

Time  tag  (GMT)  for  navigation  data 

Pseudorange  and  range  rate  (integrated 
carrier  VCO  commands)  to  satellites 
currently  being  tracked 

Navigation  data  error  variances  from 
the  Kalman  filter 

C/N0  for  all  satellites  being  tracked 

Indication  of  which  satellites  are 
being  tracked  at  any  specified  time 

Receiver  health/mode  status 

GDOP  of  primary  constellation 


Output  Buffer  2: 

GPS  satellite  almanac  data  (when 

requested) 

Output  Buffer  3: 

GMT  at  a  1  Hz  (nominal)  time  mark 
HARDWARE  DESCRIPTION 
The  GPSRU  consists  of  four  modules: 

•  RF  Module 

•  Digital  Module 

•  Interface  Module 

•  Power  Supply  Module 

A  block  diagram  of  the  GPSRU  is  shown 
in  Figure  3.  The  RF  module  includes  the 
RF/IF  translator,  the  frequency  synthe¬ 
sizer  and  the  frequency  standard.  The 
digital  module  contains  the  signal  pro¬ 
cessors,  AAMP ,  MFI  and  memory.  The 
interface  module  contains  an  AAMP,  MFI, 
memory  and  user  interface  circuitry. 

The  overall  dimensions  of  the  receiver 
shown  in  Figure  4  are  7.70  in.  wide  x 
8.28  in.  long  x  4.75  in.  high,  includ¬ 
ing  flanges.  The  unit  weight  is  8.8 
lbs.  and  the  total  power  dissipation  is 
11.6  watts.  The  size,  weight  and  power 
are  summarized  in  Table  4. 

The  modules  plug  directly  into  a  Master 
Interconnect  Board  (MIB)  for  their 
electrical  interface  and  are  secured  to 
the  chassis  with  wedgelocks.  Each  mod¬ 
ule  is  uniquely  keyed  to  prevent  cross¬ 
plugging.  A  chassis  assembly  showing 
the  plug-in  modules  is  provided  in 
Figure  5.  The  Host  Vehicle  interface 
connectors  are  mounted  on  one  side  of 
the  receiver  chassis,  as  shown  in 
Figure  4.  The  OMV  interface  consists  of 
one  9  pin  input  power  connector,  one 
SMA  type  coaxial  RF  input  connector, 
and  one  50  pin  interface  connector.  In 
addition,  there  is  one  SMA  type  coaxial 
Test  Connector. 
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PROGRAM 
MEMORY 
(80Kx  16  EPROM) 

TEMPORARY 
MEMORY 
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MEMORY 
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Host 
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Table  4.  GPSRU  Size,  Weight  and  Power  Summary 


COMPONENT 

SIZE  (in.) 

L  x  W  x  H 

WEIGHT 

(lbs) 

POWER 

DISSIPATION 

(w) 

RF  Module 

7.00  x  5.75  x  1.25 

1.10 

2.6 

Digital  Module 

7.00  x  5.75  x  0.6 

0.54 

2.6 

Interface  Module 

7.00  x  5.75  x  0.6 

0.54 

3.5 

Power  Supply  Module 

7.00  x  5.75  x  1.15 

1.25 

2.9 

MIB 

5.75  x  3.50  x  0.125 

0.34 

N/A 

RCVR  Chassis 
(Includes  Hardware) 

8.28  x  7.70  x  4.75 

5.04 

N/A 

TOTAL 

8.28  x  7.70  x  4.75 

8.8 

11.6 

SPECIFICATION 

10.25  x  8.0  x  8.5 

10.0 

18.0 

4.75 
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Figure  5.  GPSRU  Cross  Section  Side  View 


SUMMARY 

A  spaceborne  GPS  receiver  unit  (GPSRU) 
utilizing  an  advanced  technology  chip- 
set  has  been  described.  This  receiver 
weighs  8.8  pounds,  occupies  less  than 
0.2  cubic  foot  and  requires  11.6  watts. 
It  will  provide  GMT  and  spacecraft  posi¬ 
tion  and  velocity.  The  navigation  posi¬ 
tion  and  velocity  accuracy  of  the  GPSRU 
is  expected  to  be  approximately  16  feet 
(1  sigma)  per  axis  and  0.05  feet/second 
(1  sigma)  per  axis,  respectively. 
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Abstract 

This  paper  investigates  the  rise  of  GPS  as  an  aid  in  establishing  and  maintaining 
coherence  of  a  highly  distributed  Space-based  Tethered  Array  Radar  (STAR)  system.  The 
STAR  concept  is  a  network  of  distributed  small  satellite  subarrays  each  of  which  is  a  linear 
array  of  elements  or  “string”.  Each  subarray  is  made  up  of  many  dipole  elements,  each 
with  its  own  receive  module.  In  order  to  achieve  beamforming  and  the  radar  performance 
desired,  the  relative  position  and  orientation  of  each  dipole  element  must  be  estimated 
within  a  fraction  of  a  wavelength  of  the  radar  frequency. 

A  subset  of  the  general  problem  is  considered  here.  A  single  tether,  subjected  to  en¬ 
vironmental  disturbances,  is  modeled  by  a  lumped  mass  dynamic  model  of  the  motion. 
Candidate  GPS  receiver  elements  with  code  and  phase  measurements  are  located  at  the 
end- satellites  and  along  the  tether.  Based  on  this  configuration  a  linear  Kalman  filter  was 
developed  using  a  model  of  the  longitudinal  and  lateral  motions  of  the  tether.  Torsional 
effects  were  neglected  for  this  analysts. 

The  errors  caused  by  solar  pressure,  atmospheric  drag,  and  earth-shadow  transitions 
were  investigated  for  an  altitude  of  1000  kilometers.  The  estimation  error  in  the  filter 
was  determined  by  covariance  analysis,  linear  filtering  and  nonlinear  simulations- using  the 
interactive  CAE  software  Ctrl- C® /Model- (J™ . 

1.  Introduction 

Consider  two  satellites  separated  vertically  by  a  string  several  kilometers  long.  This 
is  a  very  long  antenna  that  could  be  used  for  long-wavelength  communications.  Or,  one 
could  attach  many  receivers  along  this  string  and  have  multiple  strings  next  to  each  other 
forming  a  phased  array  radar.  Such  a  space-based  radar  can  be  made  very  light-weight. 
The  launch  and  maintenance  costs  are  therefore  low  even  for  higher  altitude  orbits.  For 
a  surveillance  radar  application  it  is  necessary  that  signals  received  along  the  tether  and 
between  tethers  are  combined  coherently  to  form  a  radar  image.  It  is  therefore  neccessary 
to  know  the  relative  position  and  orientation  of  the  receive  elements  to  a  fraction  of  the 
radar  wave  length  and  the  absolute  position  to  observation  accuracy  requirements.  This 
paper  investigates  the  use  of  GPS  as  an  aid  in  establishing  and  maintaining  coherence  of 
such  a  highly  distributed  space  based  radar  system  subject  to  enivironmental  disturbances 
of  solar  pressure  and  atmospheric  drag. 

Second  Symposium  on  GPS  Applications  in  Space,  Hanscom  AFB,  Lexington,  MA, 
October  1989. 
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The  Space-Based  Tethered  Array  Radar  (STAR)  concept  under  consideration  evolved 
from  the  DOD  need  for  an  affordable,  improved  surveillance  system  [l].  The  STAR  con¬ 
cept  is  a  network  of  distributed  small  satellite  subarrays  each  of  which  is  a  linear  array 
of  elements  or  “string”  (Figure  1).  In  order  to  achieve  the  radar  performance  desired, 
the  relative  position  and  orientation  of  each  dipole  element  must  be  estimated  within  a 
fraction  of  a  wavelength  of  the  radar  frequency.  The  tether  concept,  in  addition  of  being 
lightweight,  appears  to  provide  a  basis  for  achieving  the  power  aperture  product  needed  to 
detect  the  small  targets  of  interest.  However,  the  concept  presents  a  challenge  in  terms  of 
dynamic  stability,  pointing,  and  control.  One  aspect  of  the  problem  requiring  a  solution  is 
the  development  of  an  adequate  integrated  estimator/sensor  system  to  accurately  predict 
tether  motion. 

The  approach  taken  in  this  paper  is  to  consider  a  subset  of  the  problem  and  to  analyze 
the  dynamic  behavior  of  a  single  tether  due  to  high  altitude  disturbances.  After  develop¬ 
ment  of  a  lumped  mass  dynamic  model  of  the  motion  of  a  tether  in  space,  a  candidate  GPS 
receiver  configuration  that  makes  use  of  available  GPS  code  and  phase  measurement  data 
is  discussed.  Based  on  this  configuration  and  dynamic  model,  a  linear  Kalman  filter  was 
developed  using  a  model  of  the  longitudinal  and  lateral  motions  of  the  tether.  Torsional 
effects  were  neglected  for  this  analysis. 

Section  2  describes  the  basic  tether  dynamics,  section  3  describes  candidate  GPS 
sensor  configurations,  section  4  shows  the  simulation,  and  section  5  evaluates  GPS  based 
filters. 


Figure  1:  Space-Based  Tethered  Array  Radar. 


2.  Tether  Dynamics 

The  dynamics  of  a  tether  in  orbit  around  the  earth  may  be  conveniently  analyzed 
in  terms  of  its  steady  state  and  its  transient  behavior.  The  steady  state  behavior  is  the 
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nominal  motion  of  the  tether  in  the  average  gravity  field  of  the  earth.  The  gravity  gradient 
forces  of  this  uniform  gravity  field  hold  the  tether  nominally  in  a  vertical  position  relative  to 
the  Earth’s  surface.  The  transient  behavior  is  the  flexible  motion  and  perturbation  of  the 
tether  due  to  non-uniformity  of  the  earth’s  gravity  potential,  solar  pressure,  and  interaction 
with  the  earth’s  magnetic  field.  In  adition,  thermal  heating  due  to  solar,  earthen  and  lunar 
radiation,  magnified  by  the  rapid  transition  into  and  out  of  the  Earth  shadow  will  deform 
the  string  and  introduce  further  dynamic  perturbations. 

The  dynamics  of  the  tether  may  be  modeled  as  a  continuum  or  as  interconnected 
point  masses  or  rigid  bodies.  The  continuum  model  is  well  suited  to  determine  the  steady 
state  characteristics  of  the  string  as  well  as  the  basic  eigen-frequencies.  However,  a  super¬ 
computer  would  be  needed  to  integrate  the  resulting  partial  differential  equation  derived 
from  a  force  balance  equation  written  for  the  earth’s  orbit  perturbation  environment.  For 
this  analysis  we  shall  therefore  use  discretized  models. 

Three  types  of  discretizations  are  being  used  in  modeling  the  dynamics  of  aerospace 
systems.  They  are  finite  element,  finite  difference  and  lumped  mass  models.  Finite  element 
models  are  primarily  used  for  detailed  analysis  of  nearly  linear  dynamics.  Finite  difference 
methods  are  dominant  in  fluid  and  thermodynamics.  Lumped  mass  models  are  the  easiest 
to  implement  and  are  used  for  modeling  the  overall  characteristics  of  complete  aerospace 
systems.  In  this  analysis  we  shall  model  libration  and  bowing  -  stretching  characteristics 
by  lumped  mass  models,  generally  referred  to  as  bead  models  in  the  literature  on  tether 
dynamics. 

Consider  a  string  of  point  masses  connected  by  spring/damper  combinations,  that  are 
described  in  coordinate  frames  (Orlov  frames)  that  rotate  at  the  nominal  circular  orbits 
of  these  masses  about  the  earth  (Figure  2).  The  D’Alembert  force  of  a  point  mass  is  equal 
to  the  forces  acting  on  the  point  mass. 

II 

mt~r*  =  Fg,  +  ^  FsPij  +  Fs,  +  Fs,  +  Fp{  +  Fe,  (1) 

; 

where  ()'  denotes  the  derivative  in  inertial  corrdinates,  F?t  is  the  gravity  force,  Fsp ■  is 
the  spring  force  between  the  t-th  and  j-th  mass,  Fs,  denotes  the  solar  force,  Fp%  denotes 
the  drag  force  and  F&,  is  the  electrodynamic  force. 

The  acceleration  may  also  be  described  by 

II  BB  B 

~r ,  =  y,  +  2'w  x  ~p  i +  U?  x  {Tj  X  (  R  ,  +  7\)),  (2) 

where  ~p  ,  is  the  vector  from  the  i-th  Orlov  frame  to  the  i-th  point  mass,  ()B  is  the  derivative 
relative  to  the  Orlov  frame,  w*  is  the  angular  velocity  vector  from  the  center  of  the  earth 

to  the  i-th  Orlov  frame,  and  R ,  is  the  vector  from  the  center  of  the  earth  to  the  i-th  Orlov 
frame,  where 

B 

R  i  —  —  R\k\,  R  2  =  ~[R\  +  / 1  )^2» '  ■  and  R  t  —  0  (3) 
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Coordinizing  the  vectors  in  the  Orlov  frames  gives 


Figure  2:  Bead  Model  of  Tethered  Satellite  in  Orlov  Frames. 
The  gravity  force  is, 


F  -  *mir 
*9i  ~  r2 


where  r\  —  x\  4-  y\  +  (z  -  i2*)2  and  f,  is  the  unit  vector  along  ~r\, 


r,  = 


(6) 

(7) 


giving 


(8) 
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The  spring  damper  force  of  a  tether  element  between  mass  i  and  j  is, 


FsPij  —  Cij(rO  +  Dij((  v  j  v  ,)  ■  rtj)f,j  (9) 

for  stretching  and  small  compressions,  and  is  approximately  constant  for  buckling,  where 
Cij  is  the  spring  constant,  is  the  damping  constant,  ~r* t;  is  the  vect»_.  from  mass  t  to 
mass  j,  ~v  i  is  the  velocity  of  mass  t  relative  to  the  i-th  frame,  ~v  j  is  the  velocity  of  mass 
j  relative  to  the  j-th  frame  and  /ot>  is  the  zero-force  length  of  the  tether  element  between 
mass  i  and  mass  j.  The  zero-force  length  is  a  function  of  its  temperature  and  changes  with 
the  periodicity  of  the  solar  incident  angle  and  during  Earth  shadow  transitions. 

Solar  Radiation  Forces 

An  exposed  surface  is  subjected  to  radiation  pressure  equal  to  the  difference  between 
incident  and  reflected  momentum  flux.  The  main  sources  of  electromagnetic  radiation 
pressure  in  an  Earth’s  orbit  are, 

•  Solar  illumination  (  Fs  —  1360  W  jm2  ) 

•  Solar  radiation  reflected  by  the  Earth  and  its  atmosphere  (100  -  300  W/m2  ) 

•  Radiation  emitted  by  Earth  and  its  atmosphere 

•  Solar  radiation  reflected  by  the  Moon 

The  momentum  flux,  p,  acting  on  a  surface  element  normal  to  the  direction  of  the 
radiation  is 

P  =  7  (10) 

c 

where  F  is  the  radiation  energy  flux  and  c  is  the  speed  of  light. 

The  radiation  forces  acting  on  a  surface  element  may  be  divided  into  forces  due  to 
absorbed  radiation,  forces  due  to  reflected  radiation,  and  forces  due  to  diffused  reflected 
radiation.  With  the  absorption  coefficient,  Ca,  the  reflection  coefficient,  Cr,  the  diffusion 
coefficient,  Cj,  and  Ca  4-  Cr  +  Cj  =  1,  the  total  radiation  force  for  one  radiation  source 
may  be  expressed  by, 

F  =  —p  j  (1  —  Cr)  cos  $S  +  2(Cr  cos2  cos  0)NdA,  (11) 

J  A  3 

where  0  is  the  angle  between  the  unit  surface  normal,  N,  and  the  unit  vector,  5,  towards 
the  light  source. 

The  body  surface  is  assumed  to  be  comprised  of  spherical  end-satellites,  and  flat 
surface  elements  and  cylinders  along  the  tether.  Integrating  over  different  surface  types 
we  get,  for  a  flat  surface: 


for  a  sphere  of  radius  r: 


F  -  - p4nr 2 


1  +  V 

4  +  9°d 


S, 


and  for  a  cylinder  of  height  h: 
F  =  —p2rh 


(sin  0(1  +  \cr )  +  ^ Cd)S  -  (^Crsin0  +  ^-C d)  cos  ip Z 

O  DO  O 


(13) 


(14) 


where  Z  is  along  the  cylinder  axis  and  0  is  the  angle  between  the  cylinder  axis  and  the 

A  A 

solar  direction,  0  =  arccos(Z  ■  S).  NOTE:  In  the  simulation  additional  shadowing  effects 
were  considered. 

The  total  solar  torque  may  be  computed  by  integrating  the  torques  caused  by  the 
solar  forces  acting  on  the  spacecraft.  Figure  3  shows  the  analytically  computed  torque 
due  to  solar  pressure  on  a  2000m  tether  configuration  and  the  sun  in  the  orbital  plane. 
The  total  torque  oscillates  between  ±1.4iVm  (Units  are  in  Newton-meter).  The  resulting 
steady  state  angle  oscillates  between  ±0.04  * ,  which  corresponds  to  ±1.4m  in  orbit  offset 
of  the  lower  satellite  relative  to  the  upper  satellite. 


Sun  in  Orbital  Plane:  Ca  =  1,  Cd.=  . 05.  Cr=  85 


Figure  3:  Solar  Pressure  Torque  on  Sample  Configuration. 
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Atmospheric  Forces  and  Torques 


Even  though  the  solar  pressure  forces  are  more  than  30  times  larger  than  the  forces 
due  to  atmospheric  drag,  the  drag  forces  are  significant  since  they  are  opposite  the  orbital 
velocity  and  slowly  de-orbit  the  spacecraft.  For  the  example  configuration  at  an  altitude 
of  1000/cm  this  drag  results  in  a  de-orbit  rate  of  20  m/day.  For  this  investigation,  the  1977 
Jacchia  density  model  [5]  was  used,  which  accounts  for  different  atmospheric  temperatures 
due  to  changes  in  solar  activites. 

The  force  due  to  the  impact  of  atmospheric  molecules  on  the  spacecraft  surface  can 
be  modeled  as  an  elastic  impact  without  reflection  [3].  The  differential  force,  dfAERo,  on 
a  surface  element  dA ,  with  outward  normal  N,  is  given  by, 

df\ERO  -  ~Cd^v2(N  ■  v)vdA  (15) 

where  v  =  iF/t;  is  the  unit  vector  in  the  direction  of  IF,  p  is  the  atmospheric  density  and 
Co  <  2  is  the  drag  coefficient.  The  aerodynamic  drag  for  a  surface  perpendicular  to  ~v  is 
then 

Fd  =  -Co^-A~v\~v\  (16) 

St 


The  aerodynamic  torque  on  a  spacecraft  about  its  center  of  mass  is 


where  the  integral  is  over  the  surface  of  the  tethered  spacecraft  for  which  N  •  v  >  0,  and 
~v  =  ~v  cm  +  ~u  x~p, 

Taero  =  J  CD^-  (iV  •  (~vCm  +  "w  x  7))  (lFcm  +  Zj  x  7)  x  ~pd A  (18) 

TaERO  =  J  Cd^~-{N  •  ~V  Cm)("v* cm  X  ~p)dA 

+  /  Cd^-{N  •  ("w  x  y))(lFcm  X  ~p)dA 

r  o  - 

+  J  Cd~{N  •  ~v  cm){u>  x  ~p)  x  ~pdA 
+  J  Cd^-{N  ■  (Tj  x  ~p))(~u  x  ~p)  x  ~p  dA 
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Integrating  over  the  length  of  the  tether,  the  torque  on  the  tether  about  the  center  of 
mass  pitch  axis  becomes 


Tabro,  =  j Cdy  -  ?t,ccostf<?(/23  +  /3)  +  vlcos^(il-il) 


(20) 


where  It  is  the  distance  from  the  center  of  mass  to  the  lower  S/C,  l 2  is  the  distance  from 
the  center  of  mass  to  the  string  to  the  upper  S/C,  and  6  is  the  pitch  angle.  Additional 
perturbations  due  to  atmosphere  dynamics  were  not  considered  in  this  analysis. 


3.  Candidate  Sensor  Configurations 

A  sensor  system  is  required  to  meet  the  absolute  and  relative  positioning  requirements 
of  the  radar  antenna  elements.  Knowledge  of  the  relative  positions  of  the  radar  receive 
antenna  elements  to  an  accuracy  on  the  order  of  a  fraction  of  the  radar  wavelength  is 
required  to  support  the  phase  coherency  requirements  of  the  space  based  radar.  This 
translates  to  approximately  2cm  for  L-Band  and  10cm  for  the  low  end  of  the  UHF  Band. 
The  absolute  positioning  requirement  of  the  entire  array,  with  respect  to  the  earth,  is  on 
the  order  of  10  to  100m.  A  number  of  sensor  systems  were  considered,  including: 

-  Laser  Ranging 

-  NAVSTAR  Global  Positioning  System  (GPS) 

-  Radar’s  own  signals 

-  Laser  Diode/CCD  Array 

-  Combinations  of  the  above. 

In  this  analysis  the  NAVSTAR  GPS  was  investigated,  since: 

1.  It  appears  that  GPS  can  serve  as  the  attitude  determination  sensor  in  addition  to 
meeting  both  the  absolute  and  relative  positioning  requirements  [8) 

2.  The  GPS  satellites  will  be  in  place  along  with  complete  ground  support  by  the  time 
needed  for  space  based  radar.  Thus,  only  receivers  for  the  space  based  arrays  need  to 
be  c<  nsidered.  The  GPS  solution  is  entirely  passive. 

3.  The  GPS  is  highly  robust.  At  the  altitudes  of  interest,  10  to  15  satellites  are  in  view. 
Based  on  a  solution  using  the  best  5  satellites,  the  Positional  Dilution  of  Precision 
(PDOP)  is  less  than  1.7.  A  5  satellite  solution  was  selected  as  a  trade-off  between 
solution  redundancy  and  receiver  complexity,  however,  there  are  receivers  available 
today  that  process  up  to  12  satellites  per  solution. 

4.  The  DARPA  GaAs  GPS  chip  set  has  made  available  miniature  space  qualified  hard¬ 
ware  suitable  for  space  based  radar  applications.  By  the  time  such  components  will 
be  needed  for  SBR  applications,  it  is  expected  that  multiple  sources  will  be  available. 

5.  At  least  three  programs  are  already  underway  developing  GPS  receivers  for  precise 
motion  and  attitude  sensing  in  space  applications  (  8,  9,  10,  11]. 
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6.  The  GPS  satellites  can  be  augmented  with  “pseudolites”  on  the  earth,  if  necessary, 

to  provide  a  more  robust  mode  of  absolute  position  calibration. 

Consider  now  a  configuration  with  a  5  channel  GPS  receiver  at  the  top  spacecraft 
and  multiple  differential  phase  comparators  mounted  along  the  tether,  aided  by  the  top 
receiver  for  carrier  tracking.  The  carrier  phase  of  each  satellite  can  be  tracked  to  a  few 
degrees  of  the  GPS  L-Band  carrier  which  places  the  phase  tracking  measurement  error 
at  1cm,  one  standard  deviation.  However,  the  carrier  phase  position  solution  will  be 
implemented  in  a  relative  mode  requiring  measurements  from  two  satellites,  thus  the  error 
of  the  measurement  will  be  on  the  order  of  1.4cm,  one  standard  deviation. 

For  this  analysis,  satellites  depressed  lower  than  9  degrees  below  the  local  horizon 
were  excluded  in  order  to  eliminate  those  signals  grazing  the  earth  that  would  have  large 
IONO  and  TROPO  propagation  delay  variations.  Based  on  a  unity  gain  antenna  over  this 
elevation  range  and  omnidirectional  in  azimuth,  the  received  C/A  Code  signal  strength,  in 
terms  of  carrier-to-noise  ratio  (  C/No)  available  for  tracking,  ranges  from  32  to  55 dB  Hz , 
which  is  more  than  adequate  to  maintain  track  through  the  benign  dynamics  of  the  tether 
elements.  For  comparison,  the  C/No  observed  by  GPS  receivers  in  airborne  and  terrestrial 
applications  is  in  the  range  of  25  to  30 dB  Hz. 

In  airborne  and  terrestial  applications  the  tracking  loop  bandwidths  must  be  wide 
enough  to  accommodate  the  dynamics  which  are  more  severe  than  for  the  space  based 
radar.  A  typical  airborne  GPS  receiver  might  use  tracking  loop  bandwidths  of  2  to  5 Hz 
for  the  code  loops  and  20  to  30 Hz  for  the  carrier  loop.  Because  of  the  higher  C/No  and 
lower  dynamics  environment  of  the  space  based  radar,  it  should  be  possible  to  reduce  the 
bandwidths  to  1  Hz  or  less  for  the  code  loop,  and  to  less  than  5 Hz  for  the  carrier  loop. 
The  longer  integration  times  should  result  in  sufficient  measurement  accuracy  to  support 
the  space  based  radar  position,  velocity,  and  attitude  estimation  requirements. 

4.  Simulation  Development 

The  control  and/or  estimation  problem  requires  a  linearized  model  upon  which  to 
base  the  design.  For  this  reason,  the  interactive  CAE  program  Ctrl- C® /Model- C?™  [6,7], 
was  used  for  further  analysis.  Ctrl-C^/Model-CfrM  is  a  graphical  interactive  simulation 
program.  The  components  of  the  non-linear  simulation,  masses,  springs,  and  disturbance 
forces,  are  configured  into  the  desired  model  with  a  mouse  on  the  computer  screen.  Based 
on  this  graphical  description,  Mode/- C™  generates  the  simulation  code  in  FORTRAN  or 
C  computer  languages.  A/ode/-CrM’s  greatest  benefit  for  this  analysis  is  the  capability 
to  generate  a  linearized  model  of  a  simulation  at  a  push  of  a  button  on  the  screen.  The 
linear  model  is  automatically  placed  into  the  workspace  of  the  program  to  be  available  for 
further  analysis  including  control  and  estimation  designs. 

The  nonlinear  model,  developed  within  the  Model-C ™  framework,  is  shown  in  Figure 
4.  Masses  are  connected  by  spring/dampers  and  are  forced  by  solar  pressures  that  consider 
earth  shadow  transitions.  The  sun- vector  and  solar  pressure  computations  are  performed 
in  the  “sun  block” ,  shown  in  the  upper  left  of  the  figure. 
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NON-LINEAR  MODEL 
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Figure  4:  Simulation  Model. 
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Note  that  it  is  very  easy  to  modify  and  manipulate  the  simulation.  Blocks  and  re¬ 
gions  may  be  copied  and/or  moved  around  in  the  screen.  Parameters  may  be  changed  by 
clicking  on  the  block  and  changing  a  few  numbers  in  the  opening  menu.  Outputs  may  be 
obtained  by  placing  probes  at  the  desired  points  in  the  block  diagram.  Linear  models  may 
be  obtained  by  placing  input-output  markers  and  selecting  LINEARIZE  in  the  analysis 
capabilities  menu. 

Figure  5  shows  the  probe  outputs  of  the  simulation  for  solar  pressure  forces  on  the 
upper  S/C,  Fu.solar,  the  center  bead ,  FSsolar,  and  the  lower  S/C,  FLsolar  for  a  STAR 
string  with  solar  panel  above  the  upper  S/C  and  receiver  elements  along  the  tether.  In 
order  to  analyze  the  effects  of  solar  pressure  transients,  the  simulation  was  initialized  prior 
to  the  tether  coming  out  of  the  earth  shadow.  Figure  6  shows  the  time  delay  between  the 
sun  illuminating  the  upper  and  lower  S/C. 


Tether:  h»»  10001cm,  3  — mass  Tether 


0.00  1000.  2000.  3000.  4000.  3000.  4000.  7000. 


Figure  5:  Solar  Pressure  Forces  on  the  Upper  S/C,  Center  Bead  and  Lower  S/C. 
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Figure  6:  Solar  Pressure  Forces  on  Upper  and  Lower  S/C. 

Figure  7  compares  simulations  describing  the  tether  by  three  and  five  beads.  The 
positions  of  the  center  of  the  tether  relative  to  the  upper  satellite  in  the  in-orbit  direction, 
Xt2  —  Xu,  and  in  the  vertical  direction,  Zt2  —  Zu,  do  not  change  significantly  between 
the  two  simulations  indicating  that  a  3-mass  bead  model  may  be  used  for  basic  analysis  of 
estimator  and  controller  analysis.  However,  higher  order  models  need  to  be  used  for  the 
analysis  of  the  high-frequency  bobbing  motion. 

Two  locations  along  the  orbit  were  used  to  obtain  linear  perturbation  models.  Specifi¬ 
cally,  shadow  and  sun-lit  trim  points  were  used.  Model- C?M  numerically  perturbs  dynamic 
states  and  inputs  to  generate  linear  models  of  the  form, 

*(0  =  Fx(0  +  <?*(<)  (21, 

y(t)  =  Hx(t)  +  Ju(t) 

where  x  is  the  linear  state  space,  u  is  the  input  at  the  location  of  the  input  markers  (see 
Figure  4),  and  y  is  the  linearized  output  at  the  location  of  the  output  markers. 
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Figure  7:  Comparison  of  the  Relative  Tether  Position  for  3  and  5  Mass  Simulations. 

5.  Filter  Analysis 

From  the  separation  theorem  for  linear  systems,  the  design  of  the  estimator  is  inde¬ 
pendent  of  the  design  of  the  controller.  Though  the  bowing  and  stretching  of  the  tether 
includes  non-linear  coupling,  the  separation  of  control  and  estimation  will  be  assumed  in 
order  to  get  an  upper  estimate  of  the  performance  of  the  estimator. 

This  analysis  investigates  the  capability  and  performance  of  a  Kalman  filter  to  re¬ 
construct  the  absolute  and  relative  position  and  velocity  of  receiver  elements  from  GPS 
measurements.  This  analysis  consists  of  two  parts.  These  are: 

(1)  A  covariance  analysis  is  performed  to  investigate  the  achievable  estimation  accuracies 
for  different  GPS  accuracies. 

(2)  A  linear  filter  and  its  performance  in  reconstructing  receiver  positions  and  velocities 
from  non-linear  simulation  data,  corrupted  by  receiver  noise,  is  investigated. 

Covariance  Analysis 

The  filter  was  designed  using  a  linearized  model  of  a  Model- CfrM  based  non-linear 
simulation  of  a  string, 

*(0  =  Fx(t)  +  <?iw(t) 

y(t)  =  Hx(t)  +  v(t)  {  ] 
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where  w(t)  was  chosen  to  be  a  random  noise  of  ten  percent  of  the  solar  pressure  force  plus 
a  small  contribution  stemming  from  unmodeled  dynamics,  and  v(t)  was  selected  to  be  a 
random  GPS  receiver  position  and  velocity  measurement  error, 


E(ww,T)  =  Q6{t  —  t'), 

E{vv'T)  =  R6[t  -  t '), 

E(wv>T)  =  0. 

The  GPS  configuration  model  included  a  GPS  receiver  at  the  top  satellite  providing 
absolute  position  and  velocity  based  on  code  and  phase  measurements.  In  addition,  phase 
comparator  receivers,  located  at  the  center  and  bottom  of  the  tether,  were  included  to 
provide  very  accurate  position  and  velocity  data  (relative  to  the  topmost  satellite). 

For  the  design  of  the  filter,  the  linear  model,  obtained  by  linearizing  about  the  trim 
point,  was  converted  to  a  discrete  model  at  the  sample  rate  of  the  GPS  measurement, 
giving, 

x1+i  =  Axi  +  BiWi,  E(ww J)  =  QAt  ^ 

y,-  =  Cii  + 

where  At  is  the  sample  rate.  The  sample  rate  chosen  for  this  analysis  was  a  conservative 
At  =  2  see.  if  x~  is  the  state  estimate  before  a  measurement  is  made  and  x+  is  the  state 
estimate  after  the  measurement,  we  may  define, 

M  d=f  E{x  -  x-)(x  -  x~)T,  (24) 

and  denote  the  covariance  of  the  estimation  error  after  the  measurement  update  by, 

P  d=f  E(x  -  x+)(x  -  x+)T.  (25) 

The  measurement  update  equations  for  the  Kalman  filter  are  then  [12], 

xt  =  XT  +  LdVi  ~  Hu"), 

Li  =  MiCT[CMiCT  +  R)~l, 

Pi  =  Mi  -  LiCMi 

and  the  prediction  equations  are 

x~+1  =  Axf  +  B,u„ 

Mu.l  =  APiAT  +  BxQBf, 

After  some  transient  period,  the  covariance  equations  will  converge  towards 
values,  resulting  in  the  steady  state  optimal  filter  equations,  with 

L  =  MGt(CMCt  +  R)~1, 


(26) 


(27) 

steady  state 
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(28) 


and  M  being  the  steady  state  prediction  error  covariance.  It  is  M  that  is  a  measure  of  the 
filter  performance.  The  estimated  values  of  the  state  estimation  errors  are  then  the  square 
root  of  the  diagonal  elements  of  M,  and  the  estimated  output  estimation  errors  are  the 
square  roots  of  the  diagonal  elements  of  the  output  prediction  covariance  matrix  C MCT . 

In  order  to  evaluate  the  sensor  requirements  the  achievable  accuracies  were  evaluated 
for  diffrent  absolute  and  relative  GPS  sensor  accuarcies.  Figure  9  shows  the  position 
and  velocity  estimation  errors  for  fixed  sigmas  of  the  GPS  relative  measurement  error  of 
lcm  and  .01  cm/s,  absolute  measurement  errors  ranging  from  0.2m  to  10m  and  a  2sec 
sampling  time.  The  resulting  position  estimation  errors  are  .05  -  0.55m  in  the  in-orbit 
direction,  and  .02  -  .55m  in  the  vertical  directions.  The  velocity  estimation  errors  range 
from  .0001m/s  to  .035 m/s.  These  values  are  well  within  those  required  for  the  STAR 
system.  The  depicted  lines  are  for  different  points  along  the  tether.  The  lower  values 
are  for  the  upper  end-satellite  and  the  upper  line  is  for  the  lower  end  satellite.  Figure 
10  shows  the  relative  position  and  velocity  estimation  errors  for  fixed  sigmas  of  the  GPS 
measurements  in  the  upper  S/C  of  2m  and  .02 m/s  and  relative  measurement  errors  of  the 
GPS  phase  comparison  ranging  from  .lcm  to  5cm.  The  maximum  error  of  5cm  for  relative 
position  and  .038cm/s  for  relative  velocity  are  less  than  those  permissible  for  STAR. 

Reference  (13)  shows  that  the  estimation  errors  may  be  further  reduced  by  using  a 
higher  order  models  with  an  associated  reduced  model  uncertainty.  Additional  improve¬ 
ments  may  be  achieved  by  using  a  faster  GPS  sample  rate. 


Figure  9:  Achievable  Estimation  Accuracy. 
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Figure  10:  Achievable  Relative  Estimation  Accuracy. 

Line ar  Filtering  of  Non-Line ar  Tether  Simulation  Data 

The  CAE  package  combination,  Ctrl-C®  and  Model-C?** ,  provided  a  convenient  and 
efficient  environment  in  which  to  simulate  the  non-linear  string  dynamics  and  GPS  receiver 
element  position  and  velocities  with  a  Kalman  filter  and  compare  the  estimated  outputs 
with  truth  data.  The  GPS  sensor  data  were  simulated  by  adding  uncorrelated,  Gaussian 
measurement  noise  to  the  outputs  of  the  non-linear  simulation.  The  simulation  was  ini¬ 
tialized  in  a  trim  condition  in  the  Earth  shadow  just  before  solar  forces  begin  impinging 
on  various  elements  of  the  string.  The  outputs  at  the  probe  locations  (see  Figure  4)  were 
stored  as  “truth”  data  and  corrupted  by  assumed  white,  Gaussian  measurement  noise.  To 
reduce  the  computation  time,  the  filter  equations  may  be  put  into  the  form 


*,+i  =  A  fx{  + 
Vi  =  Cx~ , 


(29) 


where 

Af  =  A  -  ALC , 
Bf  =  Cx~, 


The  combined  propagation  and  update  filter  equations  are  now  in  the  form  of  a  discrete 
simulation  and  may  be  propagated  by  a  discrete  simulator,  i.e.,  in  Cfr/-C®.  Figures  11  - 
12  show  overlays  of  the  non-linear  simulated  “truth”  data  and  estimated  values  as  well  as 
differences  between  “truth”  and  estimated  outputs.  In  the  Figures  ()u  denotes  the  upper 
end-satellite,  ( )L  denotes  the  lower  end-satellite,  ()i...  denote  locations  along  the  tether, 
()h  denotes  estimated  values,  and  6()  denote  differences  between  “truth”  and  estimated 
values. 

Figure  11  compares  absolute  position  and  velocity  estimates  for  the  upper  S/C.  The 
absolute  position  estimation  errors  are  less  than  25 cm. 

Figure  12  compares  relative  position  and  velocity  estimates  for  the  center  of  the  tether. 
The  relative  position  estimation  errors  are  less  than  3cm  at  locations  with  GPS  phase 
comparisons  and  twice  to  three  times  that  value  for  locations  without  a  relative  GPS 
measurement.  These  results  are  similar  to  those  of  the  covariance  analysis  and  indicate 
that  the  non-linear  effects  do  not  overly  burden  the  linear  filter. 

0.  Conclusions 

For  the  radar  architectures  and  tether  dynamics  considered,  it  is  shown  that  GPS 
positioning  sensors  have  sufficient  accuracy  to  maintain  radar  signal  coherence  for  the 
anatenna  elements  mounted  along  a  single  tether.  The  analysis  considered  a  tether  of 
length  2000m  at  an  altitude  of  1000/em  subject  to  solar  pressure  and  atmospheric  distur¬ 
bances,  including  earth  shadow  transitions.  A  covariance  and  Kalman  filter  analysis  was 
performed  to  estimate  the  achievable  relative  and  absolute  position  of  the  radar  receiver 
elements.  The  results  show  that  a  2m  absolute  GPS  sensor  accuracy,  a  1cm  relative  accu¬ 
racy  of  GPS  phase  comparitors  and  a  2 sec  sample  rate  suffice  to  combine  signals  along  the 
tether  with  a  Kalman  filter  to  an  accuracy  of  2cm  and  achieve  a  25cm  absolute  positioning 
accuracy  of  the  string.  The  wave  length  of  radar  frequencies  of  interest  range  from  20cm 
(L-band)  to  70cm  (UHF).  Thus  this  estimator,  given  the  GPS  error  assumptions  and  the 
considered  disturbance  environment,  can  meet  the  position  estimation  requirements  for 
coherence,  even  for  a  L-band  radar. 

Further  analysis  is  necessary  to  investigate  the  solar  heating  effects,  higher  order 
gravity  harmonics  and  unmodeled  dynamics,  particular  torsional  effects.  Non-linear  filter 
designs,  that  include  the  bowing  coupling,  may  be  required  to  get  better  estimates  for 
receiver  elements  between  sensor  locations.  Additional  analysis  is  necessary  for  coherently 
combining  signals  of  different  strings. 

Absolute  position  accuracy  requirements  are  a  function  of  mission  requirements,  radar 
performance  and  geometry.  The  shown  achievable  absolute  positioning  accuracy  of  25 cm 
seems  to  be  sufficient  for  most  applications.  However,  a  more  detailed  total  system  analysis 
must  be  performed  to  investigate  the  subsystem  error  budget  required  to  support  coherence 
among  several  tethers  in  space. 


a.  00  1000.  2000.  3000.  4000.  5000.  6000.  7000. 


Figure  11:  Absolute  Estimates  for  Top  Spacecraft. 
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ABSTRACT 

This  paper  provides  a  description  of 
the  GPS  receiver  planned  for  Installa¬ 
tion  on  the  Office  of  Naval  Research's 
(ONR)  SPINSAT  Altimeter  (SALT)  satel¬ 
lite.  The  receiver  will  use  an  ad¬ 
vanced  chipset  developed  under  con¬ 
tract  to  DARPA  that  enables  a  5-watt, 
2-pound  unit.  The  receiver  func¬ 
tional  and  hardware  description  are 
provided.  The  navigation  solution  and 
intermediate  data  available  to  the 
SALT  processor  for  telemetry  are 
identified.  The  receiver's  environ¬ 
mental  and  performance  requirements 
are  also  provided. 

INTRODUCTION 

In  August  1989,  the  Satellite  and 
Space  Electronics  Division  of  Rock¬ 
well  International  Initiated  work  on 
a  GPS  receiver  for  the  Office  of  Naval 
Research  (ONR).  The  receiver  is  sche¬ 
duled  for  delivery  on  January  31, 
1990,  and  Is  planned  to  be  Integrated 
Into  a  Commander  Class  satellite  bus 
with  an  Intended  launch  of  Spring/ 
Summer  1990.  The  primary  payload  of 
the  bus  Is  an  altimeter.  These 
efforts  are  part  of  ONR's  Special 
Purpose  Inexpensive  Satellite 
(SPINSAT)  program,  and  this  mission 
Is  referred  to  as  the  SPINSAT 
Altimeter  (SALT) . 


PURPOSE 

The  overall  SALT  mission  Is  to  deter¬ 
mine  the  viability  of  small  inexpen¬ 
sive  satellites  to  perform  operational 
missions  of  Interest  to  the  Navy.  It 
will  be  specifically  used  to  gather 
oceanographic  data.  In  order  to  per¬ 
form  this  mission,  the  position  of 
the  satellite  must  be  accurately 
known  for  combination  with  the 
altimetry  measurements.  The  proven 
doppler  beacon  system,  used  In  other 
oceanographic  missions,  will  provide 
the  primary  precision  tracking 
system.  It  will  also  provide  a 
comparative  measure  of  the  GPS 
receiver  on-orbit  performance  which 
is  proposed  to  be  utilized  In  two 
different  modes.  The  first  mode, 
discussed  more  fully  in  this  paper, 
is  the  real-time  navigation  mode. 
This  operation  is  very  similar  to 
other  aplications  of  GPS  for  navi¬ 
gating  ships  and  aircraft,  and 
provides  accurate  three-dimensional 
information  when  sufficient  GPS  satel¬ 
lites  are  in  view  of  the  receiver. 
The  second  mode  is  to  determine  highly 
accurate  positioning  in  a  post  proces¬ 
sing  mode  by  differential  operation 
with  ground  receivers.  This  mode  has 
the  potential  of  on-orbit  positioning 
to  better  than  a  meter  in  the  radial 
direction.  An  experiment  is  being 
developed  to  evaluate  this  second 
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mode  of  operation  with  the  SALT  GPS 
receiver. 

The  SALT  satellite  will  be  placed  into 
a  circular  800  kilometer  orbit  at  an 
inclination  of  108°.  The  positioning 
Information  calculated  by  the  GPS 
Receiver  and  the  measured  GPS  data 
will  be  periodically  transmitted  to 
the  ground  for  evaluation.  The  opera¬ 
tion  of  the  receiver  will  be  indepen¬ 
dent  of  the  other  onboard  systems. 
The  mission  life  projected  for  opera¬ 
tion  of  the  satellite  is  three  years. 
The  evaluation  of  using  GPS  for  this 
mission  could  lead  the  way  for  future 
scientific  and  operational  missions, 
and  reduce  the  requirements  for 
ground  tracking  sites. 


FUNCTIONAL  DESCRIPTION 

The  receiver  will  use  a  Miniature  GPS 
Receiver  (MGR)  chipset  developed  by 
Rockwell  under  contract  to  OARPaH). 
The  chips  are  shown  In  Figure  1  and 
their  physical  parameters  are  sum¬ 
marized  In  Table  1. 


The  chipset  consists  of  a  GaAs  MMIC 
RF/IF  Translator,  a  Silicon  Frequency 
Synthesizer,  a  CMOS  Signal  Processor, 
a  CMOS  Multi-Function  Interface  and  an 
Advanced  Architecture  Micro-Processor. 
This  chipset  is  configured  with  memory 
and  a  frequency  standard  into  two  mod¬ 
ules  to  make  up  the  SPINSAT  GPS  recei¬ 
ver  described  herein. 

The  receiver  has  two  channels  that 
sequentially  track  a  selected  group  of 
GPS  satellites.  The  channels  function 
Independently  so  that  data  from  two 
satellites  can  be  processed  in  paral¬ 
lel.  During  normal  tracking,  the  two 
channels  are  allocated  a  one  second 
processing  Interval  for  each  satel¬ 
lite,  i.e.,  within  one  second  all  of 
the  real-time  functions  associated 
with  acquisition,  tracking  and  naviga¬ 
tion  are  accomplished.  At  the  end  of 
the  one  second,  the  two  channels 
reacquire  and  track  the  next  two 
satellites  in  their  respective  satel¬ 
lite  assignments. 

The  two  channels  can  operate  on  either 
Lt  (1575.42  MHz)  or  l2  (1227.60  MHz). 


Figure  1.  Miniature  GPS  Receiver  (MGR)  Chip  Set 
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TABLE  1.  MGR  CHIP  SET  SUMMARY 


CHIP  TYPE 

DEVICE 

COUNt 

— 

DIMENSIONS 

(Inches) 

IMPLEMENTATION 

TECHNOLOGY 

POWER 

(mW) 

RF/IF  Translator 

420 

0.340  x  0.350 

Gallium  Arsenide 

1500 

Signal  Processor 

19,800 

0.185  x  0.220 

1.25  Micron  Bulk  CMOS 

90 

Multifunction  Interface  (MFI) 

29,000 

0.370  x  0.370 

1.6  Micron  Bulk  CMOS 

20 

Advanced  Architecture 
Microprocessor  (AAMP) 

74,000 

0.250  x  0.261 

2  Micron  Bulk  CMOS 

90 

Frequency  Synthesizer 

350 

0.165  x  0.165 

Bipolar  Silicon 

500 

Normally  L-j  Is  used  for  general 
processing  with  L2  being  used  peri¬ 
odically  to  obtain  a  separate  set  of 
measurements  to  aid  In  correcting  for 
Ionospheric  errors. 

The  receiver  can  acquire  and  track 
either  or  both  the  Clear  Access  (C/A) 
code,  or  the  Precision  Code  (P-code) 
depending  upon  performance  require¬ 
ments.  As  a  general  rule  when  esti¬ 
mated  position  Is  not  well  known  C/A 
code  Is  used  for  acquisition  and  then 
handover  occurs  to  the  P-code;  how¬ 
ever,  If  the  estimated  position  and 
time  error  Is  within  +300  P-code 
chips  (+30000  feet),  a  satellite  can 
be  acquired  with  P-code. 

A  functional  configuration  of  the 
SPINSAT  receiver  Is  shown  In  the  flow 
diagram  of  Figure  2  which  Is  designed 
to  show  the  signal/data  flow  from  the 
antenna  Input  to  the  data  output  to 
the  host  vehicle.  The  flow  diagram 
also  shows  the  functional  partition¬ 
ing  of  receiver  elements  among  RF  and 
digital  hardware,  and  software 
processing. 


HARDWARE  DESCRIPTION 

The  two-channel  SPINSAT  receiver  con¬ 
sists  of  two  modules: 

•  RF  Module 

•  Digital  Module 

A  block  diagram  of  the  receiver  Is 
shown  In  Figure  3.  The  RF  module 
Includes  the  RF/IF  translator,  the 
frequency  synthesizer  and  the  fre¬ 
quency  standard.  The  digital  module 
contains  the  signal  processors, 
Advanced  Architecture  Microprocessor 
(AAMP),  Multi-Function  Interface 
(MFI) ,  memory  and  RS232C  Interface 
circuitry. 

The  overall  dimensions  of  the  RF  mod¬ 
ule  shown  In  Figure  4  are  7.95  In.  x 
5.75  In.  x  1.25  In.  The  overall 
dimensions  of  the  digital  module 
shown  In  Figure  5  are  7.0  In.  x  5.75 
In.  x  0.6  In.  The  combined  module 
weight  and  power  requirements  of 
approximately  2  pounds  and  5  watts 
are  summarized  In  Table  2. 
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Figure  2.  SP1NSAT  Receiver  Functional  Flow  Diagram 
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Figure  3.  SPINSAT  Receiver  Block  Diagram 
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Figure  4.  RF  Module  Assembly 


Figure  5.  Digital  Module 
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TABLE  2.  WEIGHT  AND  POWER  SUMMARY 


TABLE  4.  PSEUOO  DATA  MESSAGE 


WEIGHT 

(LBS) 

WATTS 

MODULE 

+15V 

+5V 

-5V 

TOTAL 

RF 

1  .4 

B 

1 .1 

0.5 

2.1 

Digital 

0.5 

H 

- 

2.6 

TOTAL 

1  .9 

0.5 

D 

0.5 

B 

OUTPUT  DATA 

Communications  between  the  receiver 
and  the  SALT  processor  is  accom¬ 
plished  through  a  serial  data  Inter¬ 
face  Port  compatible  with  an  RS232C 
serial  port.  The  receiver  provides 
two  messages  at  a  1  Hz  rate.  These 
two  messages  consist  of  a  Navigation 
Solution  Message  (Table  3)  and  a 
Pseudo  Data  Message  (Table  4). 


TABLE  3.  NAVIGATION  SOLUTION  MESSAGE 


Mil  Of 

MUMkU 

OA1A  ITCH 
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MAT A  Tree 

or  *om 

UNITS 

1 

, 

s 

N*A* 

TlO  COOOt 
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11 
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PI  weak 

i 

r 

i 

waak 

An  iia* 

1 

FP 

i 

» 

oic 

l 

IP 

? 

• 

NiUIm 

1 

HP 

• 

• 

Velocity 

l 

fp 

A 

m/t 

Cl*ck  Hu 

1 

IP 

1 

• 

fleck  Wilt 

IP 

} 

m/% 

■icilnr  Mt 

1 

I 

t 

*a*a 

Flger*  if  Merit 

1 

I 

1 

e* no 

Petition  freer  variance* 

l 

UP 

• 

■» 
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s 
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1 
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1 
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1 
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(•/»)* 
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1 

1 
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M 

•OUSj 

m  Meat  |  tterti  U<N  alanigot 

J  PS  IWi  li  M«Mr  ef  i*c**4a  Into  An  Wart 

1  PTC  tc  *rair  *f  ceceftfi  let#  calaeaar  rear 

1  t  Integer  (14  kit*) 

•t  i  kookla  lunar  <1 

t  kit*) 

FP  »  Meets**  Point  (•  HU  aaponant. 

H  kit*  e**ti***l 
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DATA  11 CM 

PAkAMHIS 
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Naaaar 

» 

) 

t 
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1 

kl 

* 
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Channel  1  Parameter* 

1 

1 

1 
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Channel  1  Pieueo  tange 

1 

kl 

2 

B 

1 

kl 

2 
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1 

i 

1 
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1 

ki 

B 

Channel  2  k*H*  tenge  lata 

1 

kl 

m/i 

Primary  Constellation  Srit 

4 

i 

Mona 

Primary  Ca*ut«11etla*  WOP 

1 

FP 

2 

Mona 

Mav  Miter  kata 

1 

i 

1 

Mona 

Cheat* 

1 

i 

Mona 

tom 

2S 

NOUS: 

Channel  Parameter*  IncloOa:  svik 

Ate/lratk  SUta 

Carrier  Frtoueacy  (it  nr  12) 

Loot  lyu  (C/A  or  P) 

C/a*  Measurement 

m«  niter  kata  lntluCa: 

Mhir  Convergence  State* 

Channel  1  Miurmitt  Vt*4 

Channal  2  Measurement!  us*< 

•  'it  letager  (H  kit*) 

kl:  knotl*  Integer  (12  kit*) 

|  FP;  Fleeting  Paint  (k  kite  eiponant.  74  kiu  sent  Inn) 

ENVIRONMENTAL  ANO 
PERFORMANCE  REQUIREMENTS 


The  environmental  and  performance 
requirements  described  In  this  section 
relate  to  testing  of  the  receiver  at 
Rockwell  prior  to  delivery  to  the 
spacecraft  Integrator.  The  thermal 
cycle,  thermal  vacuum  and  vibration 
tests  are  summarized  In  Table  5.  The 
performance  test  parameters  and  their 
required  accuracy  are  listed  In 
Table  6. 

SUMMARY 

A  GPS  receiver  planned  as  an  experi¬ 
ment  onboard  the  ONR  SPINSAT  Altimeter 
(SALT)  satellite  scheduled  for  launch 
in  the  Spring  of  1990  has  been  des¬ 
cribed.  This  receiver  consists  of  an 
RF  and  a  digital  module  having  a 
total  weight  of  2  pounds  and  a  power 
requirement  of  5  watts.  The  unit  will 
use  an  advanced  technology  chipset 
developed  by  Rockwell  under  contract 
to  DARPA.  The  data  generated  by  the 
receiver  and  available  for  telemetry 
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TABLE  5.  RECEIVER  ENVIRONMENTAL 
TEST  REQUIREMENTS 


TEST 

REQUIREMENT 

Thermal  Cycling 

4  cycles  -  each 

(Receiver 

cycle  -1 5°C 

Operating) 

Room  Temperature 
+35°C 

Thermal  Vacuum 

2  cycles  -  each 

(Receiver 

cycle  -1 5°C 

Operating) 

+35°C 

Random  Vibration 

18G  -  RMS, 

(Receiver  Not 
Operating) 

1  minute  per  axis 

was  defined.  The  environmental  and 
performance  test  requirements  of  the 
receiver  prior  to  delivery  were 
provided. 
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TABLE  6.  RECEIVER  PERFORMANCE  TEST  REQUIREMENTS 


PARAMETER 

ACCURACY(*) 

Time-to-FIrst-FIx  at  C/Nn  =  36  dB-Hz  (LI  C/A) 

Less  than  8  min  (90%  probability) 

u 

with  Initial  PVT  uncertainty  of 

150  KM  In  position,  180  M/S  In 

velocity  and  1  sec  In  time 

Tlme-to-Achleve-NAV  Accuracy 

Less  than  2  minutes 

LI  C/A  Pseudo-Range  Error  at  C/NQ  =  36  dB-Hz 

Less  than  8  M  (1  sigma) 

LI  P  Pseudo-Range  Error  at  C/Ng  =  33  dB-Hz 

Less  than  2  M  (1  sigma) 

L2  P  Pseudo-Range  Error  at  C/Ng  =  30  dB-Hz 

Less  than  2.5  M  (1  sigma) 

LI  Delta  Range  Rate  Error  at  C/NQ  *  36  dB-Hz 

Less  than  0.03  M/S  (1  sigma) 

L2  Delta  Range  Rate  Error  at  C/Ng  *  30  dB-Hz 

Less  than  0.07  M/S  (1  sigma) 

NAV  Position  Error 

Less  than  8  M  per  axis  (1  sigma) 

NAV  Velocity  Error 

Less  than  0.15  M/S  per  axis 

(1  sigma) 

Absolute  Time  Error 

Less  than  1  microsecond  (1  sigma) 

♦Selective  Availability  Not  Enabled 
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Discussion: 

Question:  Why  Is  this  receiver  so  much  lighter  than  the  one  for  the  OMV? 

Answer:  This  receiver  does  not  have  Its  own  case  and  power  supply.  There  Is 
not  as  stringent  a  radiation  requirement  for  SPINSAT  and  shielding  is  provided 
by  the  spacecraft  housing. 

Question:  Will  the  GPS  receiver  be  run  continuously  during  the  SPINSAT  mis¬ 
sion? 

Answer:  That  will  depend  on  the  GPS  constellation  available  in  1990  and,  since 
It  Is  not  the  primary  payload,  on  whether  It  can  be  accommodated  within  the 
data  telemetry  capacity  and  other  mission  objectives. 

Question:  Will  this  receiver  be  operated  differentially  with  some  ground 
tracking  station? 

Answer:  We  are  looking  Into  that  possibility  but  It  was  not  part  of  the  origi¬ 
nal  plan. 

Question:  How  many  satellites  will  be  tracked  simultaneously? 

Answer:  The  receiver  tracks  two  simultaneously  and  switches  every  one  second 
to  maintain  track  on  a  set  of  four  primary  and  a  set  of  other  potentially 
useful  satellites.  That  makes  It  difficult  to  operate  In  the  differential 
mode  -  again,  It  was  built  to  provide  real-time  navigational  Information. 

Question:  Do  the  errors  that  you  quote  Include  satellite  segment  and  control 
segment  errors? 

Answer:  The  pseudo-range  and  range-rate  errors  are  receiver  noise  contribu¬ 
tions  only.  The  navigation  errors  Include  satellite  and  ground  control 
errors. 
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THE  DEFENSE  MAPPING  AGENCY'S  OPERATIONAL  GPS  ORBIT  PROCESSING  SYSTEM 
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DEFENSE  MAPPING  AGENCY  SYSTEMS  CENTER 


The  DMA  in  conjunction  with  NSWC  has  been  producing  global  GPS  orbits  continu¬ 
ously  for  more  than  a  decade  now  and  will  continue  to  do  so  as  funding  pro¬ 
vides.  DMA  provides  ephemerides,  tracking  data,  and  geodetic  positioning  to 
support  the  development,  test,  and  evaluation  of  weapons  systems.  By-products 
of  the  orbit  determination  include  estimates  of  satellite  clock  states  and 
earth  orientation  parameters,  as  well  as  quality  control  for  Air  Force  opera¬ 
tional  broadcast  ephemerides. 

The  GPS  tracking  network  consists  of  5  DMA  and  5  Air  Force  monitor  stations. 
These  10  permanent  tracking  stations  provide  for  global  coverage,  differential 
applications,  and  redundancy.  DMA  stations  are  outfitted  with  TI4100  recei¬ 
vers  and  a  cesium  clock;  they  collect  pseudorange,  Doppler,  navigation  mes¬ 
sage,  and  weather  data  for  daily  transmission  to  DMA.  The  receivers  are  being 
upgraded  to  handle  anti-spoofing  and  selective  availability.  Clock  analyses 
are  performed  on-site,  producing  estimates  of  frequency  bias,  effects  due  to 
aging,  and  the  Allan  variance,  In  order  to  Identify  station  clock  anomalies. 

The  data  flow  for  GPS  ephemerls  generation  begins  with  the  tracking  data  from 
the  Air  Force  and  DMA  monitor  stations  which  are  then  merged  and  fed  into 
DMA's  OMNIS  orbit  determination  software  which  produces  an  ephemerls  and  clock 
offsets.  A  reference  orbit  based  on  various  force  models  is  used  to  correct 
and  edit  the  observations  and  Is  then  combined  with  the  edited  observations  in 
the  multi -satellite  filter/ smoother.  The  data  consist  of  15-mln  smoothed 
ranges  from  the  10  stations  and  are  processed  to  produce  8-day  fits  providing 
7-day  arcs  with  half-day  overlaps  on  either  end.  Estimated  parameters  Include 
orbital  elements,  radiation  pressure  parameters,  satellite  clock  and  station 
clock  parameters,  polar  motion,  and  rate  of  change  of  UT1-UTC.  Orbital  preci¬ 
sions  are  typically  about  1  m  in  radial,  and  2  -  3  m  In  along-track  and  cross¬ 
track  components.  Along-track  error  Increases  significantly  during  eclipse 
periods.  Another  check  on  Internal  consistency  that  Is  being  Introduced  is 
the  repeated  computation  of  absolute  positions  of  each  tracking  station  using 
the  broadcast  and  precise  ephemerides.  Typical  values  In  north,  east,  and  up 
components  of  position  error  are  2  -  4  m  using  the  broadcast  ephemerides  and 
0.5  -  1.3  m  using  the  precise  ephemerides  generated  by  OMNIS. 

Future  Improvements  being  planned  Include  tuning  the  filter  for  Block  II  sat¬ 
ellites,  estimating  earth  orientation  rates,  adjusting  for  station  height 
biases,  and  possibly  accommodating  GPS  cross-link  ranging. 

Discussion: 

Question:  How  good  are  the  polar  motion  and  UT1-UTC  estimates? 

Answer:  For  polar  motion,  about  11  cm  In  x  and  y;  for  UT1-UTC,  about  0.13 
msec/day.  These  are  with  respect  to  the  IERS  final  values. 
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DMA'S  INTEREST  IN  PRECISE  GPS 

ORBITS 


*  HISTORICAL  USE  OF  TRANSIT  SATELLITES 
FOR  GEODETIC  POSITIONING  WORLDWIDE 

*  CONTINUED  NEED  FOR  GLOBAL  POSITIONING 

*  TEST  AND  EVALUATION  SUPPORT  OF 
WEAPONS  SYSTEMS 

*  SUPPORT  OF  OTHER  DMA  PRODUCTS 
(E.G.  MAPS) 
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BY-PRODUCTS  OF  DMA'S  ORBIT 
DETERMINATION 


*  PRECISE  SATELLITE  CLOCK  PARAMETER 
ESTIMATES 

*  QUALITY  CONTROL  FOR  AIR  FORCE 
BROADCAST  EPHEMERIDES,  CLOCKS  AND 
URE'S  (SPS  AND  PPS) 

*  EARTH  ORIENTATION  VALUES 
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EPHEMERIS  GENERATION  AT  DMA 


*  TRACKING  STATION  NETWORK 

*  DMA  MONITOR  STATIONS 

*  MONITOR  STATION  CLOCK  ANALYSIS  (CHAP) 

*  DATA  PREPROCESSING 

*  ORBIT  COMPUTATION  (OMNIS) 

*  ORBIT  QUALITY  ASSESSMENT 

*  ORBIT  PRODUCTS  FOR  DISTRIBUTION 
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CONTRACTOR 


APPLIED  RESEARCH  LAB. 
UNIVERSITY  OF  TEXAS 


SUPPORT  TO  DMA 


MONITOR  STATION 
HARDWARE/SOFTWARE 


NAVAL  SURFACE  WARFARE  ORBIT  DETERMINATION 

CENTER 
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GPS  TRACKING  NETWORK 

(Operational  Control  -  Segment) 

DMA  GPS  TRACKING  NETWORK 


DMA  GPS  TRACKING  NETWORK 


DMA  STATION  LOCATIONS  CHOSEN 
TO  - 


1.  ADD  HIGHER  LATITUDE  COVERAGE  IN 
NORTHERN  AND  SOUTHERN  HEMISPHERES 

2.  OBTAIN  2  (OR  MORE)  STATION 
TRACKING  OF  EACH  SATELLITE  AT  ALL 
TIMES 

*  ALLOWS  USE  OF  DIFFERENTIAL 
DATA  TYPES 

*  PROVIDES  REDUNDANCY  IN  CASE  OF 
STATION  OUTAGE 

3.  BE  ACCESSIBLE  TO  U.S.  PERSONNEL  FOR 
OPERATIONS  AND  MAINTENANCE 
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DMA  MONITOR  STATIONS 


INSTALLATION  DATES 

ARGENTINA  DEC  1985 

ENGLAND  JAN  1986 

AUSTRALIA  MAR  1986 

ECUADOR  FEB  1987 

BAHRAIN  JUN 1987 
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DMA  MONITOR  STATION  FEATURES 


*  DATA  EQUIVALENT  TO  AIR  FORCE  OPERATIONAL  CONTROL 
SEGMENT  STATIONS 

*  Tl  4100  RECEIVERS  WITH 

ARL:UT  "CORE"  SOFTWARE  IN  RECEIVER 
ARL:UT  "BEPP"  SOFTWARE  IN  PC 

*  HIGH  PRECISION  CESIUM  CLOCK 

*  COLLECT  - 

PHASE-SMOOTHED  PSEUDORANGE  (PPS) 

CONTINUOUS  COUNT  INTERPOLATED  DOPPLER  DATA 
NAVIGATION  MESSAGE 
SURFACE  WEATHER  DATA 

*  ON-SITE  QUALITY  EVALUATION 

*  DAILY  DATA  TRANSMISSION  TO  DMA 

*  ON-SITE  DATA  ARCHIVING 
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HARDWARE  BLOCK  DIAGRAM 


DUAL  GPS  RECEIVER  POWER  SUPPLY 


TI-4100  GPS  RECEIVER  #1 


TI-4100  GPS  RECEIVER  #2 


BLANK  PANEL 

HP5061A  CESIUM  FREQUENCY  STANDARD 

BLANK  PANEL 

HP  5087 A  DISTRIBUTION  AMPLIFIER 

BLANK  PANEL 


EQUIPMENT  RACK  BLOWER 


SGD 

CB890720-1 

CB 


EQUIPMENT  RACK  #  2 
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SOFTWARE 


AUTOMATED  DATA  COLLECTION  LOOP 

•  SCENARIOS  GENERATED  FROM  TEMPLATE 

•  DATA  STORED  ON  DISK 


DAILY  ACTIONS 

•  OPERATOR  INITIATED 

•  DATA  QUEUED  IN  RECEIVER  DURING  DAILY  ACTIONS 

•  QUALITY  ANALYSIS 

•  CLOCK  ANALYSIS 

•  FORMAT  DATA  FOR  TRANSMISSION 

•  MARK  DATA  FOR  ARCHIVING 

•  TRANSMIT  DATA 

•  ARCHIVE  DATA 
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r  \ 

QUALITY  FILES 


'  QUALITY  FILE  > 
TRANSMISSION 
v  PACKET 


MERGE  DATA 
INTO 

TRANSMISSION 

PACKETS 

(MERGE.BAT) 


OBSERVATION  DATA 
TRANSMISSION 
PACKET 


TRANSMIT  DATA 
TO  GE  MARK  III 
(MIRROR) 


SEND/RETRIEVE 
MESSAGE  TRAFFIC 
(MIRROR) 


GE 

MARK  III 


ARCHIVE  DATA 


FIGURE  3 


DATA 


SGD 

BT88J  220-1 
JLW 


PROCESSING  AND  TRANSMISSION 
(OPERATOR  CONTROL) 
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ON  SITE  ANALYSIS 


DATA  QUALITY  ANALYSIS 

•  SUMMARY  BY  SATELLITE:  TRACKING  STATUS 

PERCENTAGE  OF  GOOD  DATA  COLLECTED 

CLOCK  ANALYSIS 

•  INPUT  DATA:  1 5  MINUTE  AVERAGED  RANGE  RESIDUALS 

DATA  ALIGNED  IN  SIDEREAL  TIME 
RESIDUALS  BY  SATELLITE 

•  SAVE  SELECTED  DATA  IN  32  DAY  HISTORY  FILE 

•  PRODUCE  ESTIMATES  OF:  FREQUENCY  BIAS 

AGING 

STABILITY  (ALLAN  VARIANCES) 
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GPS  MONITOR  STATION  DATA  COLLECTION 

PERFORMANCE 


S 


8  8 


8  8  8  8 


o 

CM 


lN303dd  *  8861-  dSldVflO  HldflOd  dOd  031031100  VIVO 
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AUSTRALIA  ARGENTINA  ENGLAND  BAHRAIN  ECUADOR 


CLOCK  HISTORY  ANALYSIS  PROGRAM 

(CHAP) 

OBJECTIVE:  IDENTIFY  STATION  CLOCK 
PROBLEMS 

*  COMPUTES  AND  STORES  RECEIVER  TIME 

BIASES  AT  EACH  RANGE  MEASUREMENT  TIME 

*  REJECTS  BAD  DATA 

*  FIXES  LARGE  JUMPS  IN  TIME  BIAS 

*  ESTIMATES  CLOCK  PARAMETERS  FOR 

1  DAY  (CURRENT) 

2  DAYS  (CURRENT  AND  PREVIOUS) 

32  DAYS  (HISTORY) 

INCLUDING  FREQUENCY  BIAS 

AGING 

STABILITY  (ALLAN  VAR.) 


3S€» 


CHAP  CLOCK  ANALYSIS 
BASED  ON  2  DAYS 
ARGENTINA  11  JAN  89 

FREQUENCY  BIAS  BASED  ON  1-DAY  DIFFS. 


NUMBER  OF 

FREQUENCY  BIAS  (NS/DAY) 

PRN 

DIFFERENCES 

2-SIG.  AVG. 

SIGMA 

6 

13 

77.3 

2.5 

9 

22 

90.0 

4.2 

11 

2 

84.2 

1.4 

12 

19 

86.0 

2.8 

13 

17 

86.8 

4.2 

ALL 

73 

85.7 

5.4 
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CHAP  CLOCK  ANALYSIS 
BASED  ON  33  DAYS 
ARGENTINA  11  JAN  89 


LINEAR  REGRESSION  ON  TIME  BIAS  TIME  SERIES 

(UNITS:  TIME  BIAS  =  NS,  FREQ.  BIAS  =  NS/DAY) 


PRN 

TIME  BIAS 

SIGMA 

FREQ.  BIAS 

SIGMA 

3 

2835.7 

22.1 

87.7 

.7 

6 

2797.2 

11.1 

86.5 

.2 

9 

2809.9 

7.8 

86.7 

.1 

11 

2747.2 

7.3 

86.1 

.2 

12 

2800.2 

13.3 

86.7 

.2 

13 

2804.2 

10.6 

87.0 

.2 

TIME 

(DAYS) 

1 

2 

4 

8 


SORT  (ALLAN  VARIANCE) 

_ (SEC/SEC1 

2.1  E- 13 
6.9E-14 
2.6E-14 
1. IE-1 4 
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CHAP  33  DAY  CLOCK  HISTORY 


■  ■  i  ■»_  ■-  -i-  i-  ■>.  ■  ■  i  i-i 

Rmnnrini  mpnra 

*  SV  PRN  3  * 
************* 


3-DAY  WALKING  FREQUENCY  (NS/OAY) 
(SV-GROUND) 

VS.  DAY  OF  YEAR 

81.0  * 


SV-GROUND  TIME  BIAS  (NS) 

VS.  DAY  OF  YEAR 
478433- NS  &  -6638  NS/DAY  RMVD 


-620.0 


75001.0 

I 

i 


•1321.0 

1- 

348 


*********** 


******* 


***** 


************* 
*  SV  PRN  6  * 


3-DAY  WALKING  FREQUENCY  (NS/DAY) 
(SV-GROUND) 

VS.  DAY  OF  YEAR 

-1253.0 


SV-GROUND  TIME  BIAS  (NS) 

VS.  DAY  OF  YEAR 

-420000  NS  &  -1288  NS/DAY  RMVD 


174.0 


****  ** 


•1300.0 

i 

i 


-1347.0 
.  348. 


***★  *  * 


43.0  ** 

I 

i* 

i* 

i 

-88.6 
!  — 
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DMA  PRECISE  EPHEMERIS  GENERATION  DATA  FLOW 


AUSTRALIA 


ORBIT  DATA  PROCESSING 


PREPROCESSING: 

*  COMBINE  AF  AND  DMA  DATA  INTO 

IDENTICAL  FORM 

*  MODIFY  "NOMINAL  CLOCKS  AND 
EVENTS  FILE"  FOR  MASTER  CLOCK 
SWITCHES,  MOMENTUM  DUMPS  AND 
OTHER  CLOCK  EVENTS  NOT  ON  FILES 
TRANSMITTED  BY  AF  OR  DMA 
STATIONS 

*  VALIDATE  WEATHER  DATA 


OMNIS 

*  CORRECTOR  *  FILTER 

*  EDITOR  *  SMOOTHER 

*  INTEGRATOR  *  PROPAGATOR 
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OMNIS  GPS  PROCESSING  SYSTEM 


INITIAL 
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SOLVE  FOR  POLYNOMIAL  CLOCKS  AND  JUMPS 
AUTOMATIC  -  USING  CONSISTENCY  OF  ADJ.  RESIDUALS 
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TRAJECTORY  INTEGRATOR 
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IAU  80  NUTATION 

1982  DEFINITION  OF  UT1 

DMA  EARTH  ORIENTATION  PREDICTIONS 
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GENERATE  BOTH  FITTED  TRAJECTORIES  AND  SV 


PARAMETERS  ESTIMATED 
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RATE  OF  CHANGE  OF  UT1 -UTC 


DMA  GPS  ORBIT  ACCURACY 
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Along-Track  2.6  m 

Cross-Track  1.7  m 


SV06  OVERLAP  RMS  DIFF.  FOR  1988 

(In  M) 
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3  OVERLAP  RMS  DIFF.  FOR  1988 


(*)  33N3&3 JJIG 
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GPS  WEEK  n 

RADIAL  o  ALONG -TRK  A  CROSS-TRK 


SV13  OVERLAP  RMS  DIFF.  FOR  1989 

(in  M) 
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GPS  WEEK  # 
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ABSOLUTE  POSITIONING 


ALBROOK,  PANAMA 
DMA'S  GASP  SOFTWARE 

NO.  OF  INDEPENDENT  SOLUTIONS  =  16 


EPHEM./CLOCK  STD.  DEV.  OF  SOLUTION  <M) 
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DMA  GPS  EPHEMERIS  PRODUCTS 


*  PRECISE  EPHEMERIDES 

(ECEF,  POS.,  VEL.,  WGS-84) 

*  PRECISE  CLOCK  STATE  ESTIMATES 

(SATELLITE  CLOCKS) 

*  EARTH  ORIENTATION  VALUES 

*  SMOOTHED  OBSERVATIONS 

(15  MIN  CORRECTED  RANGES) 


ERS 


MEAN  S.D. 
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FUTURE  ORBIT  PROCESSING 
IMPROVEMENTS 


*  FILTER  TUNING 

BLOCK  II  SATELLITE  CLOCKS 
ORBIT  MODELING  (RAD.  PRESS.) 

*  PARTITIONING  OF  FILTER 

*  EARTH  ORIENTATION  RATES 

*  STATION  HEIGHT  ADJUSTMENT 

*  CROSS-LINK  RANGING 
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A  SHUTTLE  EXPERIMENT  TO  DEMONSTRATE  6-DEGREE-OF-FREEDOM 
NAVIGATION  WITH  GPS 

Duncan  B.  Cox,  Jr.  and,  Steven  Gardner 
MAYFLOWER  COMMUNICATIONS  COMPANY,  INC. 

and 

Neal  Carlson 
INTEGRITY  SYSTEMS,  INC. 


The  objective  of  this  effort  is  to  define  an  experiment  to  demonstrate  the 
utility  of  GPS  in  determining  position  and,  particularly,  attitude  in  space. 
The  observations  used  to  accomplish  this  are  accumulated  phase  and  pseudorange 
from  multiple  antennas.  The  unique  features  of  a  space  application  include 
high  velocity,  greater  visibility  of  GPS  satellites,  rapidly  changing  GDOP 
(many  opportunities  for  good  geometry),  and  special  environmental  effects 
(temperature,  pressure,  vibration  due  to  boost  environment).  One  goal  of  the 
experiment  is  to  achieve  low  cost  -  use  existing  equipment  and  be  able  to 
retrieve  it  for  use  after  the  experiment  and  perform  post-mission  data  analy¬ 
sis.  The  Shuttle  meets  many  of  these  requirements  for  a  low  cost  experiment 
since  it  Is  a  vehicle  specifically  designed  to  deploy  experiments  in  space. 

In  order  to  determine  attitude  one  needs  three  or  more  antennas  observing  at 
least  two  satellites,  one  In  the  plane  defined  by  the  antennas  and  one  over 
this  plane.  Although  the  Shuttle  does  have  antennas,  there  are  no  preampli¬ 
fiers;  hence  the  proposed  experiment  Includes  Its  own  deployed  array  of  anten¬ 
nas  In  the  cargo  bay.  A  determination  of  attitude  requires  an  accurate  refer¬ 
ence.  Errors  contributing  to  the  reference  Include  star-tracker  errors,  iner¬ 
tial  navigation  system  drift,  alignment  transfers,  structural  bending  between 
the  star-tracker  and  the  experiment  package.  Based  on  rotation  data  from  a 
recent  Shuttle  OMV-type  experiment,  the  difference  between  attitudes  regis¬ 
tered  by  two  Indepenent  INS's  show  a  stability  of  about  +/-  0.2  mrad  over  2000 
sec.,  which  Is  consistent  with  our  goal  of  a  few  mllllradlan  accuracy  in  the 
experiment. 

Ideally,  the  experiment  can  be  packaged  In  an  existing  experiment  module,  such 
as  the  "GAS  can"  (Hitchhiker  program)  which  is  specifically  available  for 
space  experiments.  There  are  several  options  for  mounting  the  module  In  the 
shuttle  bay  area  and  considerable  analyses  have  been  done  on  thermal  condi¬ 
tions  Inside  the  can  for  different  Shuttle  attitudes  with  respect  to  the  sun. 

The  experiment  Is  very  simple,  but  various  integration  procedures  would  delay 
actual  flight  of  the  experiment  to  about  3  years. 
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GEOMETRY 

FOR  ATTITUDE  DETERMINATION  WITH  GPS 


Mayflower  Communications 


DIFFERENTIAL  PHASE  MEASUREMENTS 
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SIX-DOF  GPS  NAVIGATION  EXPERIMENT 

Antenna  Options 
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SIX-DOF  GPS  NAVIGATION  EXPERIMENT 
Reference  Attitude  Accuracy 
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Hitchhiker-G  in  a  Typical  Flight  Configuration. 
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HITCHHIKER  G  OFFERS  USERS  THREE  MOUNTING  OPTIONS 
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Hitchhiker-M  Conceptual  Configuration. 


GET  AWAY  SPECIAL 
SMALL  SELF-CONTAINED  PAYLOADS 

TYPICAL  ORBITAL  THERMAL  ATTITUDE 
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GET  AWAY  SPECIAL 
SMALL  SELF-CONTAINED  PAYLOADS 

5  FT3  CONTAINER  WITH  INSULATED  END  CAP 


CONTAINER  TEMPERATURE  (°C» 


CUSTOt€R  PAYLOADS  PAYLOAD  PORT 
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NOMINAL  INFORMATION  RATE  OF  ONE  STANOARD  ASYNCHRONOUS  CHANNEL 


SUMMARY  OF  HITCHHIKER-G  PAYLOAD  CRITERIA 


C/5 

< 

a 

re 

CD 

To 

Cl 

05 

c 

c 

o 

o 

E 

CO 

CD 

"O 

3 

CD 

c 

c 

CO 

o 


o 

o 

Q) 

i_ 

CO 

To 

n 

T5 

C 

CO 

CO 

_o 

'c 

_o 

’> 

CO 

CD 


3 

SI 

CO 

CD 

.c 


O  CD  O 

‘-c  .E  r 

$  «3  O 

c  2  « 

CD  — 
C 

o  ~  co 

00 

E  ^  -o 

o  <  3 
o(3  w 


CO 


■D  O 


CO 

O 


—  c 

E". 

2  £ 

® « $ 

3  C/5  JO 

■3  ro  v- 

S  CO  2 
u.  a 
_  o 

P  >  -D 

TO  >  co 

■D 

I  " 

CO 


co 

E 


co 

CO 
3 
O 
O 
CO 

CO  .£ 
Q)  co 
■O  ^ 

•>  w 
o  O) 

CL  £ 


CO 

0 

To 

Q. 

■D 

05 

C 

3 

O 

E 

i 

a> 

■o 

co 


c 

o 


CO 

■o 

CO 
O 
■O  c 


CO 
TJ 

s| 

co  Q. 

0  O 
c  o 
•=  N 

CO  I 


co 

■3 

C 

3 

o 

a 

o 

Is*. 


JD 

E 

CD  co 
CO  O 
CO  •- 
<  c 
.2 
o  E 
o  a 

Q  35 

•D  ~ 
0  T3 
N  3 

2  c 


CO 


o  52 

a  0 


—  c 
o  — 

CO 

■e  -a 

05  C 

I  o 

5  a 


So 
1_  c 
0  - 
co  CO 

o  3  3 

r  2  o 

I  aa 
O  o 
rmo 

—  CM  — 


re 


E 

3 


O 
CD 

■o 
0 

O 

lO  « 
iv.  C/5 


S  ,  . 

CO  0  0 

o  _re  _re 

£  3  CL 

£  0  — 
rare 

o  re  E 

O  -I  C/5 


x 

re 


co 

0 

JZ 

cd 

£ 

LO 

05 


0 

0 

E 

re 

TO 

CO 


sz 

05 


CO 

0 

n 

o 

£ 

co 
co  <d 
CO  0£ 
CM  JC  o 
I  o  C 

lL  .e  - 


0 

c 

're 


o 

CD 


*-  X 


05 

ere 

X 


c 

o 

CD 


O  CO 

o  cx 


O  0  05 

3  re  re 

6  Q.  Q. 

CO  15  QJ  — 
C  3  05  CO 


.2  ° 

W  LO 

0 

E  ■ 


re  E 

_J  C/5 


434 


Center  of  Gravity  (CG): 

■  No  farther  than  10"  from  mounting  plate  when  accommodating  250  pounds 
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LARGE  SPACE  STRUCTURE  DISPLACEMENT  SENSING  USING  ADVANCED  GPS  TECHNOLOGY 

Gaylord  K.  Huth  and  Sergei  Udalov 
AXIOMATIX 


Aside  from  space  navigation  and  attitude  determination,  other  space  applica¬ 
tions  of  GPS  include  time  transfer  and  differential  tracking  control  of  near 
co-orbiting  spacecraft,  and  structure  control  of  large  spacecraft  (flexure 
control).  One  motivation  for  using  GPS  for  any  of  these  other  applications  Is 
reduced  cost  compared  to  competing  techniques.  This  paper  addresses  the  last 
of  the  above  mentioned  application,  l.e.  structure  control. 

Consider  a  beam  anchored  at  one  end  with  GPS  antennas  at  each  end.  Movement 
can  be  monitored  by  Interferometrlcally  observing  the  phases  of  the  GPS  carri¬ 
er  signal  at  the  two  antennas.  In  space  the  beam  Is  not  anchored,  and  the 
deflections  of  the  beam  are  referred  to  some  average  position.  The  number  of 
antennas  is  proportional  to  the  length  of  the  beam. 

A  simulation  was  done  with  a  beam  In  orbit  on  the  scale  of  the  space  station. 
The  receiver  model  assumed  the  TOPEX  receiver  and  the  multipath  contribution 
was  computed  on  basis  of  TOPEX  receiver  antennas  situated  about  0.3  m  above 
the  beam.  Multipath  was  the  major  error  contribution  In  the  simulation.  The 
beam  had  a  1  Hz  vibration  with  a  5  cm  amplitude.  The  mean  error  for  modeling 
the  vibration  using  the  differential  GPS  phase  measurements  was  about  0.6  mm, 
which  was  due  almost  entirely  to  multipath;  and  the  standard  deviation  was 
about  5  mm,  again  Influenced  mainly  by  multipath. 

It  Is  noted  from  this  simulation  that  GPS  Is  a  viable  technique  for  structure 
control,  but  there  are  other  available  techniques.  With  a  laser  retro-reflec¬ 
tive  array  or  with  Inertial  sensors  (accelerometers  and  gyros)  one  obtains 
high  performance  and  moderate  to  low  cost.  Therefore,  If  GPS  cost  can  be 
reduced  then  it  becomes  a  competitive  technique. 

Discussion: 

Question:  How  much  model  error  did  you  Introduce  Into  your  simulation? 

Answer:  The  vibration  was  assumed  to  be  perfectly  known.  We  did  not  try  to 
demonstrate  the  absolute  observability  of  actual  beam  vibrations,  but  merely 
the  capability  to  observe  a  given  vibration. 

Question:  Why  Is  GPS  being  considered  for  this  application  when  simple  accel¬ 
erometers  and  strain  gauges  will  accomplish  the  same  thing? 

Answer:  We  wanted  to  determine  If  GPS  could  do  this,  In  case  It  becomes  cost 
effective  (e.g.,  If  a  spacecraft  already  has  a  multitude  of  antennas  and  re¬ 
ceivers  on  board  for  other  purposes). 
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DISPLACEMENT  WAVEFORM 


DEFLECTION  RELATIVE 
TO  "AVERAGE"  BEAM 
POSITION 


MONTE-CARLO  SIMULATION  OVER  AN  ORBIT 
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ACTUAL  SA  EFFECTS  OBSERVED  IN  1989 

—  from  M.  S.  Braasch,  1989 
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DELAY-MULTIPLY  PSEUDORANGE  DATA  FROM  SERIES-X 
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CODELESS  GPS  TECHNIQUES 
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DELAY-MULTIPLY  RECOVERY  OF  CODE  PHASE 
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Codeless"  Delay  Processing 
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Rogue  TEC  Results:  Group  Delay  &  Phase  Delay  Comparison, 
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Rogue  TEC  Results:  Phase  Delay  Comparison, 
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